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ABSTRACT 

We present a pilot survey of 21 lines of sight towards ultracompact H II (UCH II) regions and 
three towards continuum-free lines of sight in the formaldehyde (H2CO) lio — In (6 cm) and 2n — 
2i2 (2 cm) transitions, using the H2CO centimeter lines as a molecular gas densitometer. Using 
Arecibo and Green Bank beam-matched observations, we measure the density of 51 detected H 2 CO 
line pairs and present upper limits on density for an additional 24 detected 1 10 — In lines. We analyze 
the systematic uncertainties in the H2CO densitometer, achieving H2 density measurements with 
accuracies ~ 0.1 — 0.3 dex. The densities measured are not correlated with distance, implying that it 
is possible to make accurate density measurements throughout the galaxy without a distance bias. We 
confirm that ultracompact HII regions are associated with, and possibly embedded in, gas at densities 
n(H 2 ) > 10 5 cm -3 . The densities measured in line-of-sight molecular clouds suggest that they consist 
of low volume filling factor (/ ~ 10 -2 ) gas at high (n(H2) > 10 4 cm -3 ) density, which is inconsistent 
with purely supersonic turbulence and requires high-density clumping greater than typically observed 
in gravo- turbulent simulations. We observe complex line morphologies that indicate density variations 
with velocity around UCH II regions, and we classify a subset of the UCH II molecular envelopes as 
collapsing or expanding. We compare these measurements to Bolocam Galactic Plane Survey 1.1 
mm observations, and note that most UCH II regions have 1.1 mm emission consisting of significant 
(5-70%) free-free emission and are therefore not necessarily dominated by optically thi n dust emission 
as is o ften assumed when computing clump masses. A comparison of our data with the Mang um et al. 
( |2008| ) starburst sample shows that the area filling factor of dense (n(H 2 ) ~ 10 5 cm - " 5 ) molecular 
gas m typical starburst galaxies is < 0.01, but in extreme starburst galaxies like Arp 220, is ~ 0.1, 
suggesting that Arp 220 is physically similar to an oversized UCH II region. 
Subject headings: ISM: HII regions — ISM: molecules — stars: formation 



1. INTRODUCTION 

Massive stars are known to form prefe rentially in clus- 
tered environments ( |de Wit et al.||2005 ). They therefore 
likely form from "clumps," collections of gas and dust 
more dense and compact than Giant Molecular Clouds 
(GMCs) but larger and more diffuse than typical low- 
mass protostellar cores. "Clumps" have been observed 
with masses ranging from 10 — 1O 6 M (but more typi- 
cally 10 2 — 1O 3 M ) and with beam-averaged densities in 
the range 10 3 < n(H 2 ) < 10 5 cm" 3 and sizes ^ 1 pc (e.g., 
|Rosolowsky et al.||2010| (Dunham et"aLl|2010[ ). While gi- 
ant m olecular clouds in th e Galaxy have been surveyed 
(e.g., I Jackson et al. 2006), the process by which these 
clouds condense into clumps and cores and the mecha- 
nisms by which they are dispersed are not understood. 

It is still not known what sets the final mass of massive 
stars, but it is thought that they m ust ignite while still 
accreting flMcKee fc Ostriker||2007| ). Hot O and B stars 
emitting strongly in the ultraviolet will ionize their sur- 
roundings, creating density-bounded H II regions. They 
progress from hypercompact through ultracompact and 
compact and finally diffuse H II region phases, during 
which they e ither dissociate or blow out their surround- 



ing medium ( |Churchweli'[2002] |Keto|2007[ ) . The brightest 
sources in the Galactic plane m both the free- free contin- 
uum in the cm- wavelength regime and the dust continum 
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in the sub-mm to mm-wavelength regime generally host 
UCH II regions. 

While the gas within UCH II regions is hot and ion- 
ized, the surrounding gas is initially molecular. At the 
interface between the molecular cloud and the ionization 
front, a photon-dominated or photodissociat i on re gion 
appears ( [Roshi et al.|2005[ ). |Churchwell et al.| ( [2010] ) ob- 
served HCO + towards a sample ol UCH II regions and 
noted both infall and outflow motions in molecular trac- 
ers towards these objects. It should be possible to de- 
termine whether the UCH II regions still have collapsing 
envelopes (infall signatures) or only disks (outflow signa- 
tures) and thereby determine relative evolutionary states 
of the regions. 

Two centimeter transitions of formaldehyde, 0-H2CO 
lio — In (6 cm) and 2n — 2 i2 (2 cm|M have been used 
to measure the density of molecular clouds in massive- 
star-formin g regions (e.g., |Dickel et al.| 
Goss 1987), high-latitude Galactic clo uds 
let al. 11989ft, the Galactic Center 



1986 



D ickel fc 
Turner 



(e.g., 
Zyl kt et al.|1992p 

starburst galaxies (e.g., Mangum et al.pD08|, and molec 



ular c louds in a gravitational lens (e.g., Zeiger & Darling 



20 101). Studies s i milar to o ur own have been p erformed 
"IWadi 



by | Wadiak etaT| ( |1988| ) and |Henkel et al.|fll983| ), in which 



1 All references to H2CO in this paper, except where otherwise 
noted, are to the ortho 0-H2CO population, as no para P-H2CO 
lines were observed 
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bright continuum sources were observed in the same tran- 
sitions with (approximately) beam-matched telescopes at 
~2' resolution. Our study delves deeper into the spectral 
line profiles and systematic uncertainties of H2CO den- 
sitometry and is performed at higher spatial resolution 
than past work. 

This paper presents a pilot studv as a proof-of-concept 
for a much larger ongoing surve}£] towards 400 lines of 
sight and the methodology applicable to the larger sur- 
vey. 

In section [2] we present the new observations and de- 
scribe other data sets used in our analysis. Section [3] 
describes the modeling procedure used to derive density 
from the H 2 CO line observations. Section [4] presents de- 
tailed discussion of the modeling and derivation of phys- 
ical parameters and their uncertainties. Section |5| de- 
scribes the derived and measured values. Section |6T dis- 
cusses the larger implications of our results. We conclude 
with a brief summary of important results. 

2. OBSERVATIONS AND DATA 

2.1. Source Selection 

The ob served lines-of-sight included 21 sources selected 
from the Araya et al. (2002) UCH II sample and 3 from 



the Araya et al. (2004) "massive-star forming candidate" 
sample. The sources were selected primarily on the ba- 
sis of having been previously observed with Arecibo^] 
in the lio — In transition of H2CO with the intent 
of demonstrating the densitometry method within the 
Galaxy rather than making syste matic observat i ons of 
any source class. Nonetheless, the Araya et al. (2002) 



sample includes the majority of the bright UCH II re 
gions accessible to Arecibo. Additionally, there are many 
detected GMCs along the line of sight to these UCH II 

regions. 

The I Araya et al.| ( |2004[ ) observations included 15 point- 
ings towards infrared dark cloud (IRDC) candidates 
and High-Mass Protostellar Object (HMPO) candidates. 
The sources we selected from this sample include two 
sources classified as IRDCs based on MSX data and one 
HMPO candidate. The selection of these sources was ar- 
bitrary; we were only able to observe 24 lines-of-sight in 
our 4 hour observation block. The remaining sources will 
be discussed in a later paper. The observed lines of sight 
are listed in Table Q] 



Table 1 

Measured H2CO lio — In line properties 



Source Name a 



6cm Continuum 

(Jy) 



Peak 

(Jy) 



Center 
(km s- 1 ) 



FWHM 
(km s- 1 ) 



RMS 

(Jy) 



Channel Width 
(km s- 1 ) 



G32 


80+0.19 


0.1904 


32.7968 


2 


18 


(0.01) 


-0 


393 


(0 


008) 


15.39 (0.05) 


6.57 (0.06) 


0.0049 


1 


1374 


G32 


80+0.19 1 


0.1904 


32.7968 


2 


18 


(0.01) 


-0 


092 


(0 


008) 


11.45 (0.26) 


10.25 (0.65) 


0.0049 


1 


1374 


G32 


80+0.19 2 


0.1904 


32.7968 


2 


18 


(0.01) 


-0 


063 


(0 


008) 


80.63 (0.13) 


2.49 


(0.36) 


0.0049 


1 


1374 


G32 


80+0.19 3 


0.1904 


32.7968 


2 


18 


(0.01) 


-0 


254 


(0 


008) 


84.61 (0.02) 


1.37 


(0.06) 


0.0049 


1 


1374 


G32 


80+0.19 4 


0.1904 


32.7968 


2 


18 


(0.01) 


-0 


090 


(0 


008) 


88.66 (0.09) 


3.21 


(0.31) 


0.0049 


1 


1374 


G33 


13-0.09 


-0.0949 


33.1297 





49 


(0.00) 


-0 


192 


(0 


007) 


75.92 (0.05) 


3.80 


(0.12) 


0.0045 


1 


1374 


G33 


13-0.09 1 


-0.0949 


33.1297 





49 


(0.00) 


-0 


023 


(0 


007) 


81.62 (0.35) 


2.49 


(0.88) 


0.0045 


1 


1374 


G33 


13-0.09 2 


-0.0949 


33.1297 





49 


(0.00) 


-0 


040 


(0 


007) 


101.50 (0.40) 


11.30 


(0.80) 


0.0045 


1 


1374 


G33 


13-0.09 3 


-0.0949 


33.1297 





49 


(0.00) 


-0 


039 


(0 


007) 


10.39 (0.08) 


2.04 


(0.24) 


0.0045 


1 


1374 


G33 


92+0.11 


0.1112 


33.914 





83 


(0.00) 


-0 


081 


(0 


008) 


107.28 (0.18) 


6.62 


(0.34) 


0.005 


1 


1374 


G33 


92+0.11 1 


0.1112 


33.914 





83 


(0.00) 


-0 


079 


(0 


008) 


106.03 


(0.06) 


2.41 


(0.23) 


0.005 


1 


1374 


G33 


92+0.11 2 


0.1112 


33.914 





83 


(0.00) 


-0 


160 


(0 


030) 


57.30 


(0.40) 


10.60 


(0.80) 


0.005 


1 


1374 


G34 


26+0.15 


0.1538 


34.2572 


5 


57 


(0.01) 


-1 


828 


(0 


015) 


60.24 


(0.01) 


3.80 


(0.03) 


0.0063 


1 


1374 


G34 


26+0.15 1 


0.1538 


34.2572 


5 


57 


(0.01) 


-0 


160 


(0 


015) 


26.69 


(0.08) 


1.04 


(0.22) 


0.0063 


1 


1374 


G34 


26+0.15 2 


0.1538 


34.2572 


5 


57 


(0.01) 


-0 
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(0 


015) 


11.25 


(0.19) 


2.01 


(0.40) 


0.0063 


1 


1374 


G34 


26+0.15 3 


0.1538 


34.2572 


5 


57 


(0.01) 


-0 


126 


(0 


015) 


51.70 


(2.00) 


4.20 


(1.00) 


0.0063 


1 


1374 


G34 


26+0.15 4 


0.1538 


34.2572 


5 


57 


(0.01) 


-0 


047 


(0 


015) 


48.20 


(2.00) 


1.80 


(1.00) 


0.0063 


1 


1374 


G35 


20-1.74 


-1.7409 


35.1997 


5 


17 


(0.00) 


-1 


018 


(0 


008) 


43.37 


(0.01) 


3.67 


(0.02) 


0.0051 


1 


1374 


G35 


20-1.74 1 


-1.7409 


35.1997 


5 


17 


(0.00) 


-0 


147 


(0 


008) 


36.67 


(0.10) 


1.49 


(0.27) 


0.0051 


1 


1374 


G35 


20-1.74 2 


-1.7409 


35.1997 


5 


17 


(0.00) 


-0 


324 


(0 


008) 


14.08 


(0.01) 


0.93 


(0.03) 


0.0051 


1 


1374 


G35 


20-1.74 3 


-1.7409 


35.1997 


5 


17 


(0.00) 


-0 


039 


(0 


008) 


50.59 


(0.53) 


4.92 


(1.31) 


0.0051 


1 


1374 


G35 


57-0.03 


-0.0306 


35.5779 





47 


(0.00) 


-0 


064 


(0 


009) 


52.10 


(0.10) 


4.60 


(0.30) 


0.0053 


1 


1374 


G35 


57-0.03 1 


-0.0306 
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47 


(0.00) 


-0 


021 


(0 


009) 


45.60 


(0.30) 


1.90 


(0.60) 


0.0053 


1 


1374 


G35 


57-0.03 2 


-0.0306 
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47 


(0.00) 


-0 


019 


(0 


009) 


57.60 


(0.50) 


2.90 


(0.97) 


0.0053 


1 


1374 


G35 


57-0.03 3 


-0.0306 


35.5779 





47 


(0.00) 


-0 


031 


(0 


009) 


12.80 


(0.20) 


1.84 


(0.41) 


0.0053 


1 


1374 


G35 


57-0.03 4 


-0.0306 


35.5779 





47 


(0.00) 


-0 


031 


(0 


008) 


29.04 


(0.11) 


0.82 


(0.25) 


0.0053 


1 


1374 


G35 


58+0.07 


0.0657 


35.5801 





53 


(0.01) 


-0 


146 


(0 


004) 


49.37 


(0.21) 


5.33 


(0.34) 


0.0048 


1 


1374 


G35 


58+0.07 1 


0.0657 


35.5801 





53 


(0.01) 


-0 


049 


(0 


013) 


53.13 


(0.25) 


2.98 


(0.64) 


0.0048 


1 


1374 


G35 


58+0.07 2 


0.0657 
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53 


(0.01) 


-0 


025 


(0 


004) 


58.12 


(0.29) 


3.63 


(0.74) 


0.0048 


1 


1374 


G35 


58+0.07 3 


0.0657 
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53 


(0.01) 


-0 
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(0 


004) 


13.24 


(0.17) 


2.80 


(0.39) 


0.0048 


1 


1374 


G37 


87-0.40 


-0.3993 


37.873 


4 


40 


(0.01) 


-0 


531 


(0 


006) 


60.23 


(0.11) 


8.73 


(0.35) 


0.0069 


1 


1374 


G37 


87-0.40 1 


-0.3993 


37.873 


4 


40 


(0.01) 


-0 
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(0 


014) 


53.27 


(0.19) 


4.03 


(0.46) 


0.0069 


1 


1374 


G37 


87-0.40 2 


-0.3993 


37.873 


4 


40 


(0.01) 


-0 


356 


(0 


019) 


65.13 


(0.04) 


2.74 


(0.15) 


0.0069 


1 


1374 


G37 


87-0.40 3 


-0.3993 


37.873 


4 


40 


(0.01) 


-0 


324 


(0 


045) 


72.18 


(0.04) 


1.35 


(0.14) 


0.0069 


1 


1374 


G37 


87-0.40 4 


-0.3993 


37.873 


4 


40 


(0.01) 


-0 


424 


(0 


013) 


73.97 


(0.13) 


3.01 


(0.22) 


0.0069 


1 


1374 


G37 


87-0.40 5 


-0.3993 


37.873 


4 


40 


(0.01) 


-0 


185 


(0 


012) 


79.98 


(0.06) 


1.80 


(0.14) 


0.0069 


1 


1374 


G37 


87-0.40 6 


-0.3993 


37.873 


4 


40 


(0.01) 


-0 


114 


(0 


015) 


91.96 


(0.08) 


1.21 


(0.18) 


0.0069 


1 


1374 



project code GBT10B-019 3 The Arecibo Observatory is part of the National Astronomy 

and Ionosphere Center, which is operated by Cornell University 
under a cooperative agreement with the National Science Founda- 
tion. 
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Table 1 

Measured H2CO lio — In line properties 
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a Sources are labeled by the line-of-sight followed by the number 
of the component identified, indexed from zero. The components 
do not follow a particular order, but are uniquely identifiable by 
their velocity, width, and amplitude. Table 2 

Measured H2CO 2n — 2i2 line properties 
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C^AK A ^_i_n nfi 9 
<or40.40-|-U.UO Z 


q 
o 


1 fi 
IO 


^n 

(^U 


^o^ 

UZ j 


n nnn 
u.uuu 


U.UUU j 


0.00 (0.00) 





00 


(o 


00) 


n nn/i ^ 

U.UU40 


p /i c /i7_i_n nc; n 
<or40.4 < -|-U.U0 U 


n 
u 


QQ 
OO 


? n 
(^u 


no\ 
Uz ) 


n i o/i 

-U. Iz4 


( c\ nnq^ 
(^U.UUo ) 


61.67 (0.07) 


5 


85 


(o 


17) 


n nn/i q 
u.uu4y 


r~* ak a 7 in nc; i 
<or40.4 / -|-U.UO 1 


n 
U 


QQ 
OO 


?n 

^u 


Uz ) 


n nnn 
-U.UUU 


?n nn7^ 
(^U.UU / ) 


0.00 (0.00) 





00 


(o 


00) 


n nn/i o 
u.uu4y 


P/iq ci in no n 
Lj4o.d1H-U.Uz U 


n 
U 


A 1 

41 


?n 

^u 


Uz ) 


n noo 
-U.Uzz 


( c\ nnq^ 
(^U.UUo ) 


18.50 (0.25) 


4 


39 


(o 


59) 


n nnqq 
U.UUoo 


P/| o fii j_ n H9 1 
<or4o.01-|-U.Uz 1 


n 
u 


A 1 
41 


? n 
(^u 


no\ 
Uz ) 


n nnn 
-u.uuu 


?n nnnS 
(^U.UUU ) 


0.00 (0.00) 





00 


(0 


00) 


n nnqq 
U.UUoo 


f^/IS £1 J-O H9 9 
<or4o.01-|-U.Uz Z 


n 
u 


/1 1 
41 


1 u 


no^ 

UZ j 


n nn^ 
-u.uuo 


fn nno^ 

^U.UUZ j 


52.50 (1.25) 


7 


47 


(2 


94) 


n nn*^ 

U.UUoo 


pen q9j_n fiS 
<orOU.OZ-|-U.OO U 


n 
u 


1 fi 
10 


^n 

(^U 


UZ j 


n ni i 

-U.Ull 


fn nn*^ 

^U.UUo ) 


26.21 (0.44) 


3 


10 


(1 


03) 


n nn^fi 

U.UUoO 


Pfin oo n i q n 

<orDU.OO-U.lo U 


n 
u 


OQ 

zy 


? n 
yU 


noS 
Uz j 


-U.U10 


?n nnqS 
(^U.UUo ) 


21.63 (0.21) 


2 


47 


(o 


50) 


n nnq 
U.UUo 


pci /i oj_ n no n 

<orOl.4o-|-U.Uy u 


q 
o 


A 
4Z 


? n 
(^u 


Uz ) 


n qnn 
-U.oUU 


?n c\c\a\ 

(^U.UU4 ) 


21.40 (0.02) 


2 


39 


(0 


04) 


n nnQ7 
U.UUo ( 


Pfin e/i n qs n 

<orOy.04-U.yO u 


n 
u 


0*3 
ZO 


1 u 


no^ 

UZ j 


n oon 
-u.zzu 


fn nno^ 

^U.UUZ j 


9.97 (0.03) 


5 


81 


(o 


08) 


n nn^i 

U.UUol 


P7n 9Q_i_i fin n 
<or < u.zy-|-i.ou u 


fi 




Ol 
Zl 


? n 
(^u 


no\ 
Uz ) 


n i 

-u. ioy 


?n nnqS 
(^U.UUo ) 


-23.52 (0.06) 


5 


36 


(0 


13) 


n nn/i fi 
U.UU40 


P7n 9Q_i_i fin i 
<or ( u.zy-|-i.ou 1 


fi 




Ol 
Zl 


? n 
(^u 


Uz j 


n nnn 
-u.uuu 


?n nnn^ 
(^U.UUU ) 


-0.00 (0.00) 





00 


(0 


00) 


n nn/i fi 
U.UU40 


P7n qqj_i p;q n 
<or ( u.oo-|-i.oy u 


z 


fiS 

oo 


?n 

(^U 


no^ 

UZ j 


i nsi 

-1 .UOl 


fn nn^^ 

^U.UUO j 


-21.17 (0.01) 


2 


95 


(o 


01) 

88) 


n nn^s 

U.UUoO 


JT) A Q onnf;i _i_q/i qp; n 
IriAo ZUU01-|-o4oO U 


n 
u 


nn 
uu 


(^u 


no\ 
Uz j 


-U.U10 


?n nnqS 
(^U.UUo ) 


11.51 (0.37) 


4 


14 


(o 


n nnqo 
U.UUoz 


G41. 74+0.10 


o 


28 


(0 


02) 


-0.014 


(0.002) 


14.36 


(0.34) 


3 


80 


(0 


80) 


0.0032 


G41.74+0.10 1 





28 


(o 


02) 


0.000 


(0.004) 


0.00 


(0.00) 





00 


(o 


00) 


0.0032 


G41.74+0.10 2 





28 


(o 


02) 


0.000 


(0.004) 


0.00 


(0.00) 





00 


(0 


00) 


0.0032 


G41. 74+0.10 3 





28 


(o 


02) 


0.000 


(0.004) 


0.00 


(0.00) 





00 


(0 


00) 


0.0032 


G41. 74+0.10 4 





28 


(o 


02) 


0.000 


(0.004) 


0.00 


(0.00) 





00 


(o 


00) 


0.0032 


IRDC 1923+13 





00 


(o 


02) 


0.000 


(0.000) 


0.00 


(0.00) 





00 


(o 


00) 


0.0032 


IRDC 1923+13 1 





00 


(o 


02) 


0.000 


(0.000) 


0.00 


(0.00) 





00 


(o 


00) 


0.0032 


IRDC 1923+13 2 





00 


(o 


02) 


0.000 


(0.000) 


0.00 


(0.00) 





00 


(o 


00) 


0.0032 


IRDC 1916+11 





00 


(o 


02) 


0.000 


(0.000) 


0.00 


(0.00) 





00 


(0 


00) 


0.0048 


IRDC 1916+11 1 





00 


(o 


02) 


0.000 


(0.000) 


0.00 


(0.00) 





00 


(0 


00) 


0.0048 


IRDC 1916+11 2 





00 


(o 


02) 


0.000 


(0.000) 


0.00 


(0.00) 





00 


(0 


00) 


0.0048 



a The Upper Limit Flag is 1 when the measurement indicated is 
a 3 — a upper limit on the 2n — 2\2 line depth when there is a 
corresponding lio — In line detection. 

b RMS in 1.011 km s" 1 channels. Table 3 

Distance, BGPS 1.1 mm, and other properties 



Source Name 


Distance 


Galactocentric 


KDA a 


Sl.lmm 


Source 


H 2 CO 


Scenario 13 






Distance 


Resolution 




Type 


Spectrum 






(kpc) 


(kpc) 




(Jy) 




Type 




G32.80+0.19 


12.9 


7.4 


far 


6.94 


UCHII 


red gradient 


2+3 


G32.80+0.19 1 


13.1 


7.6 


far 


6.94 


UCHII 


envelope 


2+3 


G32.80+0.19 2 


9.4 


5.1 


far 


6.94 


GMC 




2+3 


G32.80+0.19 3 


9.2 


5.0 


far 


6.94 


GMC 




2+3 


G32.80+0.19 4 


9.0 


4.9 


far 


6.94 


GMC 




2+3 


G33. 13-0.09 


9.6 


5.2 


far 


2.26 


UCHII 


red gradient 


2 


G33. 13-0.09 1 


9.3 


5.1 


far 


2.26 


GMC 


envelope 


2 


G33. 13-0.09 2 


7.1 


4.7 


tangent 


2.26 


GMC 




2 


G33. 13-0.09 3 


0.9 


7.6 


near 


2.26 


GMC 




2 



H 2 CO toward UCH lis 
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Table 3 

Distance, BGPS 1.1 mm, and other properties 



G33.92+0.11 


7.0 


4.6 


tangent 


3.86 


UCHII 


red gradient 


2 


G33. 92+0.11 1 


7.0 


4.6 


tangent 


3.86 


UCHII 


envelope 


2 


G33.92+0.11 2 


3.6 


5.8 


near 


3.86 


GMC 


- 


2 


G34.26+0.15 


3.6 


5.7 


near 


35.69 


UCHII 


red gradient 


2 


G34.26+0.15 1 


1.9 


6.9 


near 


35.69 


GMC 


- 


2 


G34. 26+0.15 2 


1.0 


7.6 


near 


35.69 


GMC 


- 


2 


G34.26+0.15 3 


3.6 


6.0 


near 


35.69 


GMC 


envelope 


2 


G34. 26+0.15 4 


3.6 


6.1 


near 


35.69 


GMC 


- 


2 


G35.20-l.74 


2.8 


6.3 


near 


- 


UCHII 


single 


4 


G35.20-l.74 1 


2.5 


6.5 


near 


- 


GMC 


- 


4 


G35.20-l.74 2 


1.1 


7.5 


near 


- 


GMC 


- 


4 


G35.20-l.74 3 


3.2 


6.1 


near 


- 


GMC 


- 


4 


G35.57-0.03 


10.3 


6.0 


far 


2.57 


UCHII 


single 


2+3 


G35.57-0.03 1 


10.7 


6.2 


far 


2.57 


GMC 


- 


2+3 


G35.57-0.03 2 


3.6 


5.9 


near 


2.57 


GMC 


- 


2+3 


G35.57-0.03 3 


1.1 


7.6 


near 


2.57 


GMC 


- 


2+3 


G35. 57-0.03 4 


2.0 


6.8 


near 


2.57 


GMC 


- 


2+3 


G35.58+0.07 


10.5 


6.1 


far 


1.44 


UCHII 


blue gradient 


2 


G35. 58+0.07 1 


10.3 


6.0 


far 


1.44 


UCHII 


- 


2 


G35. 58+0.07 2 


3.6 


5.8 


near 


1.44 


GMC 


- 


2 


G35. 58+0.07 3 


1.1 


7.5 


near 


1.44 


GMC 


- 


2 


G37.87-0.40 


9.4 


5.9 


far 


4.14 


UCHII 


blue gradient 


1 


G37. 87-0.40 1 


9.8 


6.1 


far 


4.14 


UCHII 


blue gradient 


1 


G37.87-0.40 2 


9.2 


5.7 


far 


4.14 


UCHII 


blue gradient 


1 


G37.87-0.40 3 


8.7 


5.6 


far 


4.14 


GMC 


- 


1 


G37.87-0.40 4 


8.6 


5.5 


far 


4.14 


GMC 


- 


1 


G37. 87-0.40 5 


8.1 


5.4 


far 


4.14 


GMC 


- 


1 


G37.87-0.40 6 


6.6 


5.1 


tangent 


4.14 


GMC 


- 


1 


G37. 87-0.40 7 


1.2 


7.5 


near 


4.14 


GMC 


- 


1 


G37.87-0.40 8 


1.1 


7.6 


near 


4.14 


GMC 


- 


1 


G37.87-0.40 9 


1.5 


7.2 


near 


4.14 


GMC 


- 


1 


G43.89-0.78 


8.3 


6.2 


far 


- 


UCHII 


blue gradient 


3 


G43.89-0.78 1 


8.6 


6.3 


far 


- 


GMC 


envelope 


3 


G45.07+0.13 


7.6 


6.2 


far 


4.26 


UCHII 


single 


2 


G45.07+0.13 1 


6.5 


6.0 


far 


4.26 


GMC 


- 


2 


G45.12+0.13 


7.4 


6.2 


far 


6.78 


UCHII 


other 


1 


G45.12+0.13 1 


7.4 


6.1 


far 


6.78 


UCHII 


envelope 


1 


G45.12+0.13 2 


1.9 


7.2 


near 


6.78 


GMC 


- 


1 


G45. 12+0.13 3 


7.4 


6.0 


far 


6.78 


GMC 


envelope 


1 


G45.45+0.06 


7.2 


6.1 


far 


3.71 


UCHII 


blue gradient 


2 


G45.45+0.06 1 


7.6 


6.2 


far 


3.71 


GMC 


envelope 


2 


G45.45+0.06 2 


1.9 


7.2 


near 


3.71 


GMC 


- 


2 


G45.47+0.05 


7.1 


6.1 


far 


3.34 


UCHII 


red gradient 


1+2+3 


G45.47+0.05 1 


1.9 


7.2 


near 


3.34 


GMC 


- 


1+2+3 


G48.61+0.02 


9.6 


7.5 


far 


2.20 


UCHII 


red gradient 


2+3 


G48. 61+0.02 1 


0.7 


8.0 


near 


2.20 


GMC 


- 


2+3 


G48. 61+0.02 2 


6.5 


6.4 


far 


2.20 


GMC 


- 


2+3 


G50.32+0.68 


2.1 


7.2 


near 


- 


UCHII 


- 


1 


G60.88-0.13 


2.8 


7.4 


near 


4.90 


UCHII 


limit 


2 


G61.48+0.09 


5.2 


7.5 


far 


7.86 


UCHII 


single 


4 


G69.54-0.98 


2.57 


7.9 


tangent 


- 


UCHII 


thick 


4+5 


G70.29+1.60 


7.3 


9.1 


far 


- 


UCHII 


blue gradient 


2 


G70.29+1.60 1 


7.8 


9.3 


far 


- 


GMC 


envelope 


2 


G70.33+1.59 


7.3 


9.1 


far 


- 


UCHII 


single 


1+2 


IRAS 20051+3435 


2.6 


7.6 


tangent 


- 


GMC 


limit 




G41.74+0.10 


11.3 


7.6 


far 


0.56 


UCHII 


limit 




G41. 74+0.10 1 


11.6 


7.7 


far 


0.56 


UCHII 






G41. 74+0.10 2 


2.4 


6.8 


near 


0.56 


GMC 






G41. 74+0.10 3 


3.8 


6.1 


near 


0.56 


GMC 






G41. 74+0.10 4 


11.2 


7.4 


far 


0.56 


UCHII 






IRDC 1923+13 


4.2 


6.5 


near 




GMC 


limit 




IRDC 1923+13 1 


5.5 


6.3 


tangent 




GMC 






IRDC 1923+13 2 


3.8 


6.6 


near 




GMC 






IRDC 1916+11 


2.0 


7.2 


near 




GMC 


limit 




IRDC 1916+11 1 


4.2 


6.2 


near 




GMC 






IRDC 1916+11 2 


3.6 


6.4 


near 




GMC 







a The Kinematic Distance Ambiguity described in Section |5.3| 
b Scenario or scenarios most likely to be consistent with the ob- 
served spectrum, as described in Section |6.3| In some cases, the 
spectrum was consistent with multiple scenarios or some blend of 
multiple scenarios. In others, the source could not be classified, in 
which case it is marked with -1 in this column. 
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Table 4 

Inferred H2CO line properties 



Source Ncime 


'l — 1 


n 1 (TFC) 


'2 — 2 


to (FFC) 
• 2 — 2 v y 


2-2 Upper 

J_jlllll I r Icxg, 


2cm Are£i a 

//2 


6cm Aresi a 
112 


FFC Error 


G32.80+0.19 


0.18 (0.055) 


0.2 (0.059) 


0.12 (0.024) 


0.15 (0.031) 





88.0 


226.2 


0.1 


G32.80+0.19 1 


0.04 (0.013) 


0.043 (0.013) 


0.016 (0.0051) 


0.021 (0.0065) 





88.0 


226.2 


0.1 


G32.80+0.19 2 


0.027 (0.0089) 


0.029 (0.0095) 


0.0033 (0.00069) 


0.0042 (0.00088) 





88.0 


226.2 


0.1 


G32.80+0.19 3 


0.11 (0.035) 


0.12 (0.037) 


0.014 (0.0028) 


0.018 (0.0035) 





88.0 


226.2 


0.1 


G32.80+0.19 4 


0.039 (0.012) 


0.042 (0.013) 


0.0055 (0.0011) 


0.0071 (0.0014) 





88.0 


226.2 


0.1 


G33. 13-0.09 


0.34 (0.1) 


0.49 (0.15) 


0.16 (0.032) 


0.63 (0.12) 





33.5 


33.5 


0.2 


G33. 13-0.09 1 


0.035 (0.015) 


0.047 (0.02) 


(0.0059) 


(0.0031) 


1 


33.5 


33.5 


0.2 


G33. 13-0.09 2 


0.062 (0.022) 


0.084 (0.029) 


(0.0059) 


(0.0031) 


1 


33.5 


33.5 


0.2 


G33. 13-0.09 3 


0.061 (0.021) 


0.082 (0.028) 


(0.0059) 


(0.0031) 


1 


33.5 


33.5 


0.2 


G33.92+0.11 


0.084 (0.027) 


0.1 (0.031) 


0.045 (0.0091) 


0.094 (0.018) 





214.0 


214.0 


0.2 


G33.92+0.11 1 


0.082 (0.026) 


0.098 (0.031) 


0.036 (0.0072) 


0.075 (0.014) 





214.0 


214.0 


0.2 


G33.92+0.11 2 


0.17 (0.062) 


0.21 (0.074) 


(0.0049) 


(0.0031) 


1 


214.0 


214.0 


0.2 


G34.26+0.15 


0.38 (0.12) 


0.4 (0.12) 


0.22 (0.043) 


0.26 (0.052) 





10.9 


10.9 


0.2 


G34.26+0.15 1 


0.028 (0.0089) 


0.029 (0.0092) 


0.0067 (0.0014) 


0.0079 (0.0017) 





10.9 


10.9 


0.2 


G34.26+0.15 2 


0.017 (0.0059) 


0.018 (0.006) 


0.0026 (0.00059) 


0.0031 (0.0007) 





10.9 


10.9 


0.2 


G34.26+0.15 3 


0.022 (0.0072) 


0.023 (0.0074) 


0.0036 (0.00092) 


0.0043 (0.0011) 





10.9 


10.9 


0.2 


G34.26+0.15 4 


0.0082 (0.0036) 


0.0085 (0.0037) 


0.0026 (0.0011) 


0.003 (0.0013) 





10.9 


10.9 


0.2 


G35.20-l.74 


0.21 (0.063) 


0.22 (0.066) 


0.071 (0.014) 


0.084 (0.017) 





39.5 


39.5 


0.2 


G35.20-l.74 1 


0.028 (0.0085) 


0.029 (0.0088) 


0.0039 (0.0011) 


0.0046 (0.0013) 





39.5 


39.5 


0.2 


G35.20-l.74 2 


0.063 (0.019) 


0.065 (0.019) 


0.008 (0.0017) 


0.0095 (0.0019) 





39.5 


39.5 


0.2 


G35.20-l.74 3 


0.0073 (0.0027) 


0.0075 (0.0028) 


(0.0023) 


(0.0031) 


1 


39.5 


39.5 


0.2 


G35.57-0.03 


0.11 (0.035) 


0.15 (0.049) 


0.056 (0.011) 


0.26 (0.054) 





6.7 


6.7 


0.1 


G35.57-0.03 1 


0.034 (0.018) 


0.046 (0.024) 


0.011 (0.0047) 


0.047 (0.02) 





6.7 


6.7 


0.1 


G35.57-0.03 2 


0.03 (0.017) 


0.042 (0.023) 


(0.0099) 


(0.019) 


1 


6.7 


6.7 


0.1 


G35.57-0.03 3 


0.05 (0.021) 


0.069 (0.029) 


(0.0099) 


(0.019) 


1 


6.7 


6.7 


0.1 


G35.57-0.03 4 


0.051 (0.02) 


0.069 (0.028) 


0.017 (0.0065) 


0.077 (0.029) 





6.7 


6.7 


0.1 


G35.58+0.07 


0.24 (0.071) 


0.32 (0.097) 


0.085 (0.017) 


0.61 (0.12) 





2.1 


2.1 


0.2 


G35.58+0.07 1 


0.072 (0.029) 


0.096 (0.038) 


(0.0072) 


(0.019) 


1 


2.1 


2.1 


0.2 


G35.58+0.07 2 


0.037 (0.012) 


0.049 (0.016) 


(0.0072) 


(0.019) 


1 


2.1 


2.1 


0.2 


G35.58+0.07 3 


0.05 (0.016) 


0.066 (0.021) 


(0.0072) 


(0.01) 


1 


2.1 


2.1 


0.2 


G37.87-0.40 


0.12 (0.037) 


0.13 (0.038) 


0.047 (0.0095) 


0.061 (0.012) 





27.5 


170.9 


0.2 


G37.87-0.40 1 


0.028 (0.0089) 


0.029 (0.0092) 


0.0095 (0.0024) 


0.012 (0.0031) 





27.5 


170.9 


0.2 


G37.87-0.40 2 


0.081 (0.025) 


0.084 (0.026) 


0.0075 (0.0021) 


0.0096 (0.0027) 





27.5 


170.9 


0.2 


G37.87-0.40 3 


0.074 (0.024) 


0.076 (0.025) 


0.011 (0.0023) 


0.014 (0.003) 





27.5 


170.9 


0.2 


G37.87-0.40 4 


0.097 (0.029) 


0.1 (0.03) 


0.0098 (0.002) 


0.013 (0.0026) 





27.5 


170.9 


0.2 


G37.87-0.40 5 


0.041 (0.013) 


0.043 (0.013) 


0.0033 (0.00068) 


0.0043 (0.00087) 





27.5 


170.9 


0.2 


G37.87-0.40 6 


0.025 (0.0083) 


0.026 (0.0086) 


0.0021 (0.00052) 


0.0026 (0.00066) 





27.5 


170.9 


0.2 


G37.87-0.40 7 


0.039 (0.012) 


0.04 (0.012) 


0.0054 (0.0012) 


0.0069 (0.0015) 





27.5 


170.9 


0.2 


G37.87-0.40 8 


0.016 (0.0069) 


0.016 (0.0071) 


0.0035 (0.00081) 


0.0046 (0.001) 





27.5 


170.9 


0.2 


G37.87-0.40 9 


0.03 (0.0095) 


0.031 (0.0098) 


0.0035 (0.00073) 


0.0045 (0.00094) 





27.5 


170.9 


0.2 


G43.89-0.78 


0.25 (0.074) 


0.32 (0.096) 


0.037 (0.0078) 


0.12 (0.024) 





13.5 


13.5 


0.1 


G43.89-0.78 1 


0.025 (0.0077) 


0.031 (0.0097) 


0.0097 (0.0022) 


0.029 (0.0067) 





13.5 


13.5 


0.1 


G45.07+0.13 


0.092 (0.029) 


0.13 (0.04) 


0.04 (0.0081) 


0.096 (0.02) 





2.5 


2.5 


0.2 


G45.07+0.13 1 


0.058 (0.02) 


0.08 (0.027) 


0.0061 (0.0019) 


0.014 (0.0045) 





2.5 


2.5 


0.2 


G45. 12+0.13 


0.043 (0.013) 


0.045 (0.013) 


0.014 (0.0028) 


0.017 (0.0033) 





15.4 


516.6 


0.2 


G45.12+0.13 1 


0.035 (0.011) 


0.036 (0.011) 


0.0095 (0.002) 


0.011 (0.0025) 





15.4 


516.6 


0.2 


G45. 12+0.13 2 


0.046 (0.014) 


0.048 (0.014) 


0.0075 (0.0015) 


0.009 (0.0019) 





15.4 


516.6 


0.2 


G45. 12+0.13 3 


0.006 (0.002) 


0.0062 (0.0021) 


0.0033 (0.00068) 


0.004 (0.00082) 





15.4 


516.6 


0.2 


G45.45+0.06 


0.32 (0.096) 


0.32 (0.095) 


0.063 (0.013) 


0.069 (0.012) 





1963.0 


1963.0 


0.2 


G45.45+0.06 1 


0.025 (0.011) 


0.026 (0.011) 


0.01 (0.0021) 


0.011 (0.0019) 





1963.0 


1963.0 


0.2 


G45.45+0.06 2 


0.011 (0.0036) 


0.012 (0.0035) 


(0.0031) 


(0.01) 


1 


1963.0 


1963.0 


0.2 


G45.47+0.05 


0.35 (0.11) 


0.45 (0.14) 


0.089 (0.018) 


0.39 (0.079) 





3.0 


3.0 


0.2 


G45.47+0.05 1 


0.018 (0.0068) 


0.023 (0.0084) 


(0.01) 


(0.01) 


1 


3.0 


3.0 


0.2 


G48.61+0.02 


0.058 (0.018) 


0.068 (0.021) 


0.015 (0.0034) 


0.053 (0.012) 





25.5 


25.5 


0.2 


G48.61+0.02 1 


0.02 (0.0075) 


0.023 (0.0088) 


(0.0067) 


(0.0026) 


1 


25.5 


25.5 


0.2 


G48.61+0.02 2 


0.016 (0.0052) 


0.018 (0.0061) 


0.0033 (0.0013) 


0.012 (0.0046) 





25.5 


25.5 


0.2 


G50.32+0.68 


0.027 (0.01) 


0.045 (0.017) 


0.0089 (0.0031) 


0.058 (0.019) 





108.0 


108.0 


0.2 


G60.88-0.13 


0.12 (0.037) 


0.14 (0.043) 


0.011 (0.0031) 


0.031 (0.0071) 





615.0 


615.0 


0.2 


G61.48+0.09 


0.088 (0.026) 


0.09 (0.027) 


0.069 (0.014) 


0.088 (0.017) 





355.0 


355.0 


0.2 


G69. 54-0.98 


0.98 (0.29) 


5.7 (1.7) 


0.18 (0.037) 


2.9 (0.57) 





0.5 


0.5 


0.2 


G70.29+1.60 


0.086 (0.026) 


0.089 (0.027) 


0.022 (0.0044) 


0.026 (0.0052) 





52.8 


52.8 


0.1 


G70.29+1.60 1 


0.011 (0.0037) 


0.012 (0.0038) 


(0.0019) 


(0.0026) 


1 


52.8 


52.8 


0.1 


G70.33+1.59 


0.7 (0.21) 


0.78 (0.24) 


0.34 (0.068) 


0.52 (0.1) 





16.4 


16.4 


0.1 


IRAS 20051+3435 


0.12 (0.036) 


0.13 (0.014) 


0.015 (0.0041) 


0.016 (0.0034) 





2747.75 


2747.75 


0.0 


G41. 74+0.10 


0.13 (0.04) 


0.2 (0.06) 


0.01 (0.0027) 


0.045 (0.012) 





75.2 


75.2 


0.2 


G41. 74+0.10 1 


0.04 (0.014) 


0.06 (0.021) 


(0.0071) 


(0.0026) 




75.2 


75.2 


0.2 


G41. 74+0.10 2 


0.14 (0.043) 


0.21 (0.065) 


(0.0071) 


(0.0026) 




75.2 


75.2 


0.2 


G41. 74+0.10 3 


0.045 (0.017) 


0.067 (0.025) 


(0.0071) 


(0.0026) 




75.2 


75.2 


0.2 


G41. 74+0.10 4 


0.089 (0.029) 


0.13 (0.043) 


(0.0071) 


(0.0028) 




75.2 


75.2 


0.2 


IRDC 1923+13 


0.02 (0.0062) 


0.02 (0.0047) 


(0.009) 


(0.0028) 




2747.75 


2747.75 


0.0 


IRDC 1923+13 1 


0.016 (0.0051) 


0.017 (0.0038) 


(0.009) 


(0.0028) 




2747.75 


2747.75 


0.0 


IRDC 1923+13 2 


0.0081 (0.0028) 


0.0083 (0.0023) 


(0.009) 


(0.0028) 




2747.75 


2747.75 


0.0 



H 2 CO toward UCH lis 
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Table 4 

Inferred H2CO line properties 

IRDC 1916+11 0.033 (0.011) 0.036 (0.0062) (0.014) (0.042) 

IRDC 1916+11 1 0.082 (0.025) 0.089 (0.0095) (0.014) (0.042) 

IRDC 1916+11 2 0.017 (0.0064) 0.018 (0.0046) (0.014) (0.042) 



1 


2747.75 


2747.75 


0.0 


1 


2747.75 


2747.75 


0.0 


1 


2747.75 


2747.75 


0.0 



a The beam area is 2747. 75 //2 , which is used when the CMB is 
the only background continuum illumination 



Table 5 

Derived physical properties from H2CO 



Source Name 


N(H 2 CO) a 


N(H 2 CO) (FFC) b 


n(H 2 ) a 


n(H 2 ) (FFC) b 


X H 2 CO a 


X H2 co (FFC) b 


Flag c 




(cm ) 


(cm J 


— 3\ 

(cm J 


(cm J 




G32.80+0.19 




12.941°]! 


5.10l°-g 


1- +0.27 
°-^ ± — 29 


-10-79to20 


-10.75l°i| 


2 


G32.80+0.19 1 


12.05+°}? 

-0.11 


12.14+° ]% 

— 0.13 


4.96 + °^ 

— 0.28 


5.05+°" 

-0.28 


-11.39+0-22 

— 0.23 


-11.39+°" 

0.20 


2 


poo on 1 n 1 O O 

Kj6Z. oU+U. 19 z 


11.66_ 10 


11 71 +0.10 
11 71 -0 .10 


z, 1fi +0.39 
4 - 16 -0.38 


z, oq+0.31 
4t.^O_ 32 


in Q7+0- 44 
-10.97_ 46 


111 n+° 37 

-ll.10_o.3T 



z 


G32.80+0.19 3 


I2.l8±g;jg 


12.23l°°9 


4-07±g;l| 


4.23l°|l 


— 10 Q7+0- 44 
lU.J/_ 45 


-10.48lg;|S 


2 


G32.80+0.19 4 


11.82+g-;° 
— u.uy 


11.87+?, 12 
— u.uy 


4.30+giJ 


4.44+°l| 

— U. 


-10.97+HI 

— U.O 1 


-11.05+21' 


2 


G33. 13-0.09 


> 12.80 


> 13.56 


> 4.54 


> 5.10 


> -10.62 


> -11.70 


8 


G33. 13-0.09 1 


< 11.96 


< 11.90 


< 4.50 


< 3.91 


< -8.44 


< -8.45 


6 


G33. 13-0.09 2 


< 12.20 


< 0.00 


< 4.29 


< 0.00 


< -8.29 


< 0.00 


5 


G33. 13-0.09 3 


< 12.20 


< 0.00 


< 4.32 


< 0.00 


< -8.29 


< 0.00 


5 


G33.92+0.11 


> 12.35 


> 12.64 


> 4.86 


> 5.16 


> -11.29 


> -12.30 


8 


G33.92+0.11 1 


12 34+ - 07 


12.65_S;S 


4 97+0.22 

' -0.23 


5.26^ 


-11 11+ - 19 

±± - ±± -0.22 


-11 09 +0 13 

- 1 - - LUtl _0 .16 


2 


G33.92+0.11 2 














9 


G34.26+0.15 


_ ^ _ _i_n 1 n 

i3.oi±g:i? 


13.13tg;|S 


. _i_n ok 
4 Q1 ^ u - zo 
^• yi -0.29 


- ~ — i— ^1 

OUi -0.32 


-10.38tg;|l 


- ~ ^ ^ — 1— fl 1 7 

-10.36lg;_ 


2 


G34.26+0.15 1 


1 1 7q^ u - uy 
— 0.08 




^• D ' -0.25 


, in 21 

4 - ' °-0.24 


—11 Qfi^ 


— 1 1 40^ 


2 


G34.26+0.15 2 


- ^-i-D 1 n 

ll.53l°;}° 


11.561°]° 


4.38l°j° 


^ ^o_0.30 


-ll.33l°;|° 


, ^-t-fl ^2 

— 1 1 4fK 


2 


G34.26+0.15 3 


ll.63l°; 11 


iiOO -0.10 


4 -^_0 .32 




-11.28l°;|| 


- - ^ , _i_n ^1 
iit)4 -0.32 


2 


G34.26+0.15 4 


- - a ^-l-D 1 7 

11.401°,; ;i 


1 1 4"^ ' 

±±4AO -0.16 


4.87lg;2J 


, ^ , If) Qfl 

4,y4 -0.42 


-ll.95l°;|° 


iJ - yo -0.34 


2 


G35.20-l.74 


12.601°;°? 


- ^. „ _ _i_n ns 
12.65l°°? 


4 - ' z -0.25 


, _^ in occ 
^ ' ^-0.26 


_ _ _i_n 9^ 
-10.61Io.27 


- ^ ^, . un 22 

-io.82lS:12 


2 


G35.20-l.74 1 


ll.69l°;}J 


libido 


4 ^n +u '°° 


4 41 tu.ou 
^ •^■ L -0.33 


— 11 ncp 
ii.uy_ 39 


- - - 1— O "\A 

— 11 1 R ' 

iiiO -0.36 


2 


G35.20-l.74 2 


1 1 Q7^ u - uy 
±x.y f _ .o9 


_ „ „ _i_n n q 

12.001°°° 


4 9fl + 
4 - ZU -0.35 


^ -1-0 *^0 
4 ^1 ^ u ou 
^•«5^_0.30 


-10.70±g;2? 


-10.79lo.te 


2 


G35.zU-l.74 3 


< 11.41 


< 11.44 


< 4.89 


< 5.02 


< —9.24 


< —9.30 


6 


G35.57-0.03 


> 12.42 


> 13.38 


> 4.82 


> 5.61 


> -11.20 


> -12.02 


8 


G35.57-0.03 1 


> 11.72 


> 12.25 


> 4.51 


> 5.13 


> -11.80 


> -12.72 


8 


G35.57-0.03 2 


< 11.98 


< 12.12 


< 4.96 


< 5.12 


< -8.71 


< -8.72 


6 


G35.57-0.03 3 


< 12.09 


< 12.23 


< 4.60 


< 4.78 


< -8.38 


< -8.37 


6 


G35.57-0.03 4 


> 11.93 


> 12.47 


> 4.58 


> 5.20 


> -11.55 


> -12.48 


8 


cv^ ^s+n 07 


> 12.58 


> 14.06 


> 4.50 


> 5.48 


> —10.79 


> —11.71 


§ 


G35.58+0.07 1 


< 12.19 


< 12.32 


< 4.08 


< 4.55 


< -8.07 


< -8.14 


6 


G35.58+0.07 2 


< 11.96 


< 12.10 


< 4.57 


< 4.95 


< -8.53 


< -8.63 


6 


G35.58+0.07 3 


< 12.06 


< 12.17 


< 4.35 


< 4.38 


< -8.32 


< -8.24 


6 


G37.87-0.40 


-1 C\ A A +0 07 

12 - 44 Io.07 


t n >~ <-»-|-0 07 

1253 l0.09 


o „-|-0 22 
4 - 8 6lo.23 


4.98^24 


-1 rx /-»,"k + 20 

-10.90^;^ 


-10.92_J-_; 


2 


oqt 07 n An 1 
Cj37.o7-U.4U 1 


1 1 07+0. 10 
-L-L-O'_0.09 


1 1 QC+0 10 

ll.b»D_ og 


A qo+0.22 
4 - 89 -0.26 


c nn +0.21 
5.00_ 24 


11 cn+O- 21 
-11.50_ 24 


1 rn+0.18 

iiOt) -0.22 


2 


G37.87-0.40 2 


12.101°,]* 


12 09 +0 14 


16 +i.i5 

0.1D_ 1 20 


7Q+1.78 
' y -0.71 


-9.5411:5? 


-10.181°;^ 


4 


G37.87-0.40 3 


12.05t°:i° 


12 in +0 10 


4 00+0.32 

^•°°-0.33 


4 49+0 .28 

^•^-0.29 


— 10 7«+°- 37 
iu. m_ 39 


-10.87lg;|i 


2 


G37.87-0.40 4 


12.171°;^ 




3-36lJ;i 


4 13+0.35 

^ -L*J_0.36 


-9 67+ 1 - 12 

y-O' -1.60 


-10.46l° ;« 


2 


G37.87-0.40 5 


11.851°;^ 


ll.88lS;_S 


3.05+i'°* 


Q 04+I 33 


-9.68±i:i§ 


-00411." 


4 


G37.87-0.40 6 


11 67+ 019 
ll.D(_ 23 


ii.rolS;_S 




qq + 1.32 

• «JO -j^ -j^ 3 


-9.89±i;|g 


— 1 n 1 1 + 1 26 

iUii -1.58 


4 


G37.87-0.40 7 


ii-8ii°o:o° 9 




4.28l°'l 




-10.95l°;g 


-11 07 +0 32 

J-J-.U ( _ 32 


2 


G37.87-0.40 8 


11-5612:11 


11.6S1S;1J 


4.67l°; 3 3 ° 


4.8otS;s; 


-11.591CJI 


-11 fip;+ 030 

AAOD -0.34 


2 


G37.87-0.40 9 


11 72+ - 13 

— 0.29 


11^ .0° 9 


3.98lJ;S 


4.32lS;g 


-10.74l° 2 .i° 


—1 1 nfi +0 36 


2 


G43.89-0.78 


124q +o.io 


12.76^;°o 7 


A ,0+0.34 

4.18_o. 33 


4.68l° ; 2 2 8 8 


-10.17l°; 3 4 J 


-10 40 +0 24 


2 


G43.89-0.78 1 




12.80lS; 6 00 


4 95 +0 20 


6 16+ 096 


-11.66lg;l» 


-11.84lg-|| 


1 


G45.07+0.13 


12.38^ 


12 7K+0.13 
' —0.20 


4.96lg;g 


c OK+0.25 
DZD -0.27 


-11 06+ 019 
ll.UD_ 22 


-io.97^; 1 1 5 8 


2 


G45.07+0.13 1 














9 


G45.12+0.13 


12 02 +0 - 08 
iz.uz_ 07 


12 H7+ 08 
ijiU '-0 .07 


4.83+° ; 2 2 l 


4 92 +0 19 

^ yZ -0.21 


-11 so + °- 2 ° 

ll.JU_ 23 


-ii.32tS:E 


2 


G45.12+0.13 1 


11 qn+0-09 
ii.yu_ .o 8 


11 qc+008 


4 74+0.22 

4 - '^-0.23 


4 00+O.2I 
4.»«5_ 23 


1 1 Qo+0.23 
ii.jz_ 25 


-ll.S6lS;S 


2 


G45.12+0.13 2 


11 90+ - 10 


11 qq+009 
iiyti -0 .09 


4 41+0- 26 

4 - 4i -0.28 




-ll.00±g;ll 


-ii06l° ; 29 


2 


G45.12+0.13 3 


11 /L£+0- 08 
ll-4«_ 09 


ii,OO -0.12 


K 1K + 0.19 

°- iO -0.21 


°^-0.22 


-12 i^+°- 16 
iz.io_ lg 


-12.16is.17 


2 


G45.45+0.06 


12.62l°;°l 


12.64lg;„ 


4.33l°; 2 3 ? 


4 Q 7 +o.28 

4t -« 3 ' -0.28 


-10.19l°|3 


-in 21 +0 30 


2 


G45.45+0.06 1 


11.89l°;°l 


11.02lg;g? 


5.001°;^ 


5 04+° 26 


-ii.59±g : i| 


—1 1 fin +0 23 


2 
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Table 5 

Derived physical properties from H2CO 



G45.45+0.06 2 
G45.47+0.05 
G45.47+0.05 1 
G48.61+0.02 
G48.61+0.02 1 
G48.61+0.02 2 
G50.32+0.68 
G60.88-0.13 
G61.48+0.09 
G69.54-0.98 
G70.29+1.60 
G70.29+1.60 1 
G70.33+1.59 
IRAS 20051+3435 
G41. 74+0.10 
G41. 74+0.10 1 
G41. 74+0.10 2 
G41. 74+0.10 3 
G41. 74+0.10 4 
IRDC 1923+13 
IRDC 1923+13 1 
IRDC 1923+13 2 
IRDC 1916+11 
IRDC 1916+11 1 
IRDC 1916+11 2 



< 11.55 

19 71 +0.09 
1Z. I 1_ 07 

< 11.91 

12 06 +0 - 09 

< 11.88 

> 11.71 

12 24 +0 - 18 

-0.25 

> 12.51 



19 91 +0.09 
iZ,Zi -0.08 

< 11.53 

i3.i6i°-;s 

12 - 20 5# 

12.25_ 23 

< 12.12 

< 12.18 

< 12.17 

< 12.27 

< 11.86 

< 11.86 

< 13.29 

< 12.05 

< 12.36 

< 12.09 



< 11.66 
13-48±S;S 

< 11.92 

12.53^1? 

< 11.91 
11 94+° 18 

' — . 23 

> 12.41 

-,0 0^+0.09 

> 12.62 



12 ok+0.08 
iZ ^°-0.08 

< 11.55 

13 41+° 19 

±0 ^ -0.35 

19 9 o+0.04 

12.48±S;SS 

< 0.00 

< 0.00 

< 0.00 

< 0.00 

< 11.84 

< 11.77 

< 11.53 

< 12.57 

< 12.40 

< 13.77 



< 4.77 

4 46+ - 28 

^•^°-0.28 

< 5.34 

4.68t°;£ 

< 5.09 
4 60 +0 - 28 

> 4.61 
3-20±i;ig 

> 5.07 



4.67^ 



7 +0.23 
-0.23 

< 4.50 

4 64 +0 - 34 

^• D ^-0.32 

2 99 +0 - 99 

< 4.72 

< 3.21 

< 4.70 

< 4.11 

< 5.19 

< 5.30 

< 8.00 

< 5.16 

< 4.61 

< 5.55 



< 5.39 



5.21 



+0.40 



-0.34 

< 5.23 

c oq+0.22 
DZb, -0.24 

< 4.47 

^ 09+0.23 

OZZ -0.36 

> 5.31 
4 51 +0 27 

> 5.19 



4 74+0 23 

^- '^-0.24 

< 4.67 
4.83tS:»? 

4.iii o: 2 23 

4 50 +0 28 

< 0.00 

< 0.00 

< 0.00 

< 0.00 

< 4.47 

< 4.58 

< 5.00 

< 5.62 

< 5.03 

< 8.00 



< -9.04 

-10.23±°;i; 

< -8.65 

-11 io +0 - 25 

n.iu_ 26 

< -8.61 

-11 VL+O- 33 

— 0.33 

> -11.77 

_q 44+ 1 - 29 

^•^-1.43 

> -11.27 



< 

-9.96^ 



-10 Q4+ - 24 

-8.92 
-+0.22 

" -0.31 

-10 40 +0 45 

y,zo -l .22 

-8.37 
-8.91 
-8.32 
-8.22 
-8.63 
-8.70 
-8.95 
-8.44 
-8.13 
-8.68 



< 
< 
< 
< 
< 
< 
< 
< 
< 
< 



< -9.43 
-10.2llS;S 

< -8.57 

-11 2^ +0 13 

— 0.16 

< -8.58 
-11 76 +0 14 

— 0.19 

> -12.63 
-10.64l° ; 2 ? 

> -12.33 

-1 07+° 23 
Au y ' -0.26 

< -8.98 
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5 
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a The values used in this paper are shown in boldface. Uncorrected 
values are listed in this column. The filling-factor corrected values 
are shown for comparison in the next column even though they 
were not used for analysis. 

k The values used in this paper are shown in boldface. Filling- 
factor corrected values are listed in this column. The uncorrected 
values are shown for comparison in the previous column even 
though they not used for analysis. 
c Flags: 

1. No filling factor correction (no FFC) is the most reliable. 

2. Filling factor correction (FFC) is the most reliable 

3. There is an ambiguity between low density / high abun- 
dance and low abundance / high density (no FFC) 

4. There is an ambiguity between low density / high abun- 
dance and low abundance / high density (FFC) 

5. Upper Limit (No FFC) 

6. Upper Limit (FFC) 

7. Lower Limit (No FFC) 

8. Lower Limit (FFC) 

9. Unreliable estimate because of continuum / filling factor 
uncertainty. 

10. No limit (S/N) 

11. Optically Thick 

2.2. Green Bank Telescope 

We observed the H2CO 2n — 2 12 line at 2 cm 
(14.488789 GHz) with the Green Bank Telescope (GBTf] 
dual-beam Ku-band receiver as part of project GBT09(> 
049. The GBT dual-beam Ku-band receiver was used for 
4 hours on January 18th, 2010 in beam- switched nodding 
mode. System temperatures ranged from 27 to 38 K in 

4 The National Radio Astronomy Observatory operates the 
GBT and VLA and is a facility of the National Science Foundation 
operated under cooperative agreement by Associated Universities, 
Inc. 



the H2CO band centered on the 2n — 2 12 line. A band- 
width of 12.5 MHz (258.8 km s _1 ) and channel width 
of 3.052 kHz (0.063 km s _1 ) were used with 9- level sam- 
pling, with receiver temperature « 2 IK. Three additional 
tunings were acquired simultaneously, centered between 
the H and He 75a, 76a, and 77a radio recombination 
lines (RRLs) with the same channel widths and band- 
widths as above at 14.1315, 14.6930, and 15.2846 GHz. 
Each source was observed for 150 seconds in each re- 
ceiver for a total on-source integration time of 300 sec- 
onds. Each observation in the pair was independently 
inspected to search for emission/absorption in the off po- 
sition, which was 5.5' away in azimuth. When absorption 
was detected in one of the off positions, that on/off pair 
was discarded if one of the detected lines was affected, 
but otherwise was noted and ignored. Pointing and fo- 
cus observations on the calibrator source 1822-0938 were 
taken at the start of and two hours into the observations. 

The gain was assumed to be 1.91 K/Jy based on pre- 
vious calibration observations on point sources in Ku- 
band; our flux density measurements will therefore be 
overestimates for extended sources. The aperture effi- 
ciency was t]a = 0.671, and the main beam efficiency 
t]m b = 1.3277,4 = 0.886, so our main-beam corrected 
measurements could overestimate extended source flux 
densities by at most 13% (ignoring atmospheric ab- 
sorption). The data were calibrated using the normal 
GETNOD procedure in GBTIDL^], which assumes an at- 
mospheric opacity at zenith r = 0.014 at 14.488 GHz. 

We assume primary beam 0fwhm = 51.1" per the 
GBT observers manual. We assume a conservative 10% 
error in the beam area Q = 7.8 x 10 -8 sr, which governs 

5 GBTIDL Jhttp://gbtidl .nrao.edu/ ) is the data reduction pack- 
age produced by JNKAO and written in the IDL language for the 
reduction of GBT data. The National Radio Astronomy Observa- 
tory is a facility of the National Science Foundation operated under 
cooperative agreement by Associated Universities, Inc. 
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the flux density received from the CMB over the observed 
area. Beam size error should be dominated by small 
errors in focus. By utilizing the 305 m Arecibo telescope 
at 6 cm and the 100 m GBT at 2 cm, we acquired beam- 
matched (FWHM- 50") observations of the H 2 CO li - 
ln and 2n — 2 12 lines. 

2.2.1. GBT Data Reduction 

In the 24 lines of sight, 75 independent components 
were identified from the lio — In spectra. These were fit 
with gaussians using GBTIDL's fitgauss routine. Out 
of these 75 components, 51 had corresponding 2n — 2i2 
detections. The fitted gaussian spectral lines are listed 
by line of sight in Table [I] The gaussian fits may not be 
representative of the true spectral line profi le; c omplex 
spectral line profiles are discussed in Section \6?2\ 

The 2 cm continua were measured by fitting a first- 
order baseline in each reduced nodded pair excluding the 
line and the bandpass edges. Figure [I] shows the flat 
baselines achieved in the observations, though the RRL 
spectrum shows an example of the artifacts seen at the 
edges of the bandpass. The continuum error listed in the 
table is the RMS of only the data included in the baseline 
fit after the baseline was subtracted from the spectrum; 
the systematic error from flux calibration uncertainty is 
20% and dominant. 

2.3. Arecibo 

The Arecibo 4.829660 GHz H 2 CO li - In observa- 
tions used in this project we re previo u sly pr esented in 
Araya et al. (2002) and Araya et al. (2004) and were 
kindly provided in reduced form by E. Araya. All obser- 
vations were performed using standard on/off position 
switching and 5 minute integration times in both the on 
and off positions, resulting in off positions 1.25 degrees 
away from the pointing center. We assume a 30% error 
in the continuum (ba sed on measured gains in the range 
2.0-2.5 as reported in | Araya et al.||2002[ ) and an effective 
diameter of 227m (Ofwhm = 56'", Vt 9.0 x 10~ 8 sr ) 
with 10% uncertainty] 

The Arecibo spectral lines were re-fit for this paper by 
converting the Arecibo data from CLAS^jto GBTIDL's 
SDFITS format^ and using GBTIDL's FITGAU SS rou- 
tine. The 6 cm continua were taken directly from |Araya| 
eTaD ([20021) Table 3. 



G32.8Q + 0.19 



2.4. Other Archival Data 
2.4.1. Very Large Array 

We acquired VLA archival images from the Multi- 
Array Galactic Pl ane Imaging Survey (MAGPIS) 6 cm 
Epoch 3 data set flHelfand et*aL||20Q6| ) and the NRAO 
VLA Archive Survey (NVAS^ MAGPIS has a resolu- 
tion of ~ 4" and sensitivity a ~ 2.5mJy/bm. The NVAS 
has variable resolution and sensitivity since it is based on 
VLA archival data. The VLA data was used to estimate 
source sizes and interferometer-to-single-dish flux ratios. 

6 http : //www . naic . edu/$\sim$phil/sysperf /misc/hpbw_vs_ 
lambda_2004.html 

' CLASS is part of the GILDAS software developed by IRAM. 

8 Code for t he CLASS-GBTIDL conversion is available from 
http : / / code . google . com/p/ casaradio/wiki/ class_to_gbt 

y The NVAS is run by Lourant Sjouwerman at the NRAO. It 
has not yet been published. 
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Figure 1. Top: The GBT 2n - 2 12 (red) and Arecibo lio - In 
(black) spectra of G32.80+0.19. Bottom: The GBT H75a (red) 
and Arecibo HI 10a (black) spectra with the GRS 13 CO spectrum 
(light blue) overlaid. The left axis is for the RRLs and the right 
axis is for the 13 CO. The C and He RRLs are not displayed. 

2.4.2. Bolocam 1.1 mm 

We extract 1.1 mm dust continuum fluxes from the 
Bolocam Galactic Plane Survey (BGPS) vl.O release 
summing over a 25" radius aperture after subtracting the 
median in a 50-200" annulus to remove background con- 
tributions. The aperture size is selected to match the 1.1 
mm data to the 2 and 6 cm data. We assume a uniform 
50% systematic error in BGPS fluxes from combined un- 
certai nties in the ca l ibrati on and background subtrac- 
tion. Aguirre et al. (2011) contains a complete discus- 
sion of the uncerta inties in the BGPS. We apply the 
Aguirre et al. (2011) recommended flux correction of 1.5 
and aperture correction for a 25" aperture of 1.21. Ad- 
ditionally, d ata from the Bolocam catalog (Rosolowsky 



et al.||2010[ ) was used with the flux correction and an 
aperture correction of 1.46 for 20" apertures. 

2.4.3. Boston University / Five College Radio Observatory 
Galactic Ring Survey 

The BU FCRAO GRS flJackson et al.|2006[ ) is a survey 
of the Galactic plane in the l6 C0 1-0 line with ~ 46" 
resolution. We extracted spectra in 25" radius apertures 
from the publicly available data for comparison with the 
H 2 CO spectra. 

2.4.4. GLIMPSE 

The Galact ic Legacy Infrared M id-Plane Survey Ex- 
traordinaire (Benjamin et al.lpOS GLIMPSE) maps 
were used to examine the morphology of the objects in 
our survey in order to determine whether an IRDC was 
present. 

3. MODELS AND ERROR ESTIMATION 

A grid of large velocity gradient (LVG) models w as run 
using both the RADEX ([van der Tak et al I|20p7| ) code 
and a proprietary code by Henkel et al. Q198(j| ) with a 
gradient of 1 km s _1 pc _1 . The models from the two 
codes were consistent to within ~ 10% in predicted op- 
tical depth and T ex . Both utilized collision rates from 
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Green 
and mil 



1991) extracted from the LAMDA databas^ 
tiplied by the recommended factor of 1.6 to ac- 



count for collisions with H2 being more efficient than 
He. The expected accuracy is ~ 30%. Zeiger fc Dar- 



|hng| |2010j) demonstrated that the errors in collision rates 
lead to systematic errors < 50% (0.3 dex) in the mea- 
sured quantities (N(H2CO), n(H2)). When measuring 
density and column, we used the RADEX models be- 
cause of their extensively tested code and documenta- 
tion. All of the models used a kinetic temperature of 40 
K and covered a range of 500 densities x 500 columns 
logarithmically sampled over 10 1 < n(H.2) < 10 7 cm -3 
and 10 11 < 7V(o-H 2 CO) < 10 16 cm" 2 . The assumption 
Tk = 40 K is reasonable in UCH II regions, which should 
be warmer than IRDCs and other cold molecular clouds. 
Dust tempera tures measured towards UCH II regions are 
around 40 K ( |Rivera-Ingraham et al.|2010| ). In the fore- 
ground clouds, this assumption is less well supported, 
but as long as the temperatures are higher than ~ 20 K, 
the models change little with temperature (Figure [2]). 

Because of a collisional selection effect, above its 
critical density (n cr (H2CO lio — In) ~ 8 cm~ 3 , 
n cr (H 2 CO 2u - 2 i2 ) « 76 cm" 3 , |Mangum et"al]|2QQ8 ) 
H 2 CO preferentially overpopulates lower st ates of the K- 
doub let (AJ 



Henkel et al. 



0, AK a = 0, AK C = ±1 
1980[ ) . These spectral lines are cooled to excitation tem- 
peratures lower than the CMB and can therefore be seen 
in absorption against it. The lio — ln/2n — 2i2 ab- 
sorption line ratio is sensitive to the density of H2 at 
densities > 10 3,5 cm -3 , allowing measurements of the 
density to within ~0.3 de x with little sens itivity to gas 
kinetic temperature (|Mangum et al.| [2008). When den- 
sity is 'measured' with critical density based tracers such 
as CO, CS, HCN, or HCO+, the estimate can be off by 
as much as 2 orders of magnitude because of radiative 
trapping effects. Similarly, measurements of density as- 
suming spherical symmetry can be very far from the local 
values. 

The collision rates of H2CO with H2 ha ve been re- 
derived with a claimed accuracy of 10 % by |Troscompt| 



et al.| ( |2009b[ ). |Troscompt et al.|(|2009a[ ) showed that col- 
lisions with para-H2 are more efficient at cooling H2CO 
into absorption against the CMB than He or ortho-H2 
and that H2CO absorption is therefore sensitive to the 
Ortho/Para ratio of H2. These improved rates are not 
used in this paper since they are only computed over a 
more limited range of temperatures, but may be used in 
future works. 

3.1. Turbulence 
Molecular gas is often observed to have spect ral line 



widths consistent with supersonic turbulence (Kainu- 



lainen et al.|2009|) and therefo re a lognormal density dis 
tribution (|Kritsuk et al.|2007| ). Our LVG models assume 
constant density per velocity bin, so the resulting mod- 
els should be smoothed by the probability distribution 
function (PDF) of the density. For clouds with a nar- 
row density distribution (logarithmic standard deviation 



of the density cr s = <J\ 



n(p)/ln{p) 



< 



smoothing is smaller than other syst< 



the effect of 
ic errors, but 



10 Ihttp : //www . strw . leidenuniv . nl/$\sim$moldata/ 

11 We use p to indicate number density in this section in order 
to be consistent with the cited literature. Because the widths are 
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Figure 2. The predicted optical depth ratio {top) and optical 
depth {bottom) vs. volume density assuming a fixed abundance 
I -H 2 CO — 10 -9 per km s _1 pc _1 shows that the dependence of 
the derived density on temperature is weak. At lower abundances, 
these curves shift to the right, providing sensitivity to moderately 
higher densities. Our 5-cr detection limit is generally around r ~ 
0.01. 

for more turbulent clouds the density PDF width can 
excee d an order of magnitude (e.g., |Federrath et al. 
2010) and will substantially change the derived density. 
Because the Mach numbers of the turbulence in the ob- 
served clouds are unconstrained, we cannot correct for 
this added uncertainty. The change in measured den- 
sity is |Alog(p)| < 0.25 for a s < 0.5, with a slight 
bias towards higher densities at lower optical depth ratio 
r iio-iii/ r 2ii-2i 2 (Figure [3]). However, for a s = 1.5, the 
bias exceeds an order of magnitude at some densities. 

Additionally, we consider the effects of "gravoturbu- 
lence", in which a high-density tail inconsistent with a 
lognor mal distribution is observed. Kainulai nen et a"L] 
(2009) report column density distributions derived from 
2 MASS extinction measurements that can be used as a 
proxy for the density distribution for a wide variety of 
clouds. Non-star- forming clouds retain a lognormal dis- 
tribution and are consistent with the analysis presented 
above. However, evolved star- forming regions develop 
a high-column density tail. For evolved (actively star- 
forming) regions like Ophiucus, Orion, and Perseus, the 
high-column density tail is substantial, and H2CO den- 
sity measurements will be highly biased towards the high- 
est density gas. More quiescent regions like the Pipe and 
Coalsack nebulae are consistent with a lognormal col- 
umn distribution to a degree that the high-column den- 
sity tail would not affect H2CO density measurements 
significantly. 

To demonstrate the effects of turbulent distributions, 
we calculate the optical depth ratio as a function of the 
mean density for three turbulent widths in Figure [3] We 
compare the density that would be inferred from the 
spectral line ratio assuming no turbulence (just LVG) to 
the 'correct' density including turbulent effects in Figure 
[4] We have also compared the LVG and turbulent densi- 
ties to "gravoturbulent" density distributions, in which a 
power law tail of high-densit y gas begins at about 10~ 2 
times the peak density (e.g., |Klessen| 2000 ; Cho & Kim 



2011 ), but because the density distributions m these sim- 
ulations are relatively narrow, the effects of the high- 
density tail on the measured density are negligible except 

relative to a mean density, the scaling between mass and number 
density is unimportant. 
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Figure 3. The optical depth ratio as a function of density for 
turbulent density distributions with widths specified in the legend. 
The optical depth ratio varies more slowly with density than in the 
pure LVG model (the solid line is the same as the black 10 K line 
in Figure |2ji). 

for the most turbulent cases. 

Figure [3] is meant to demonstrate the effects of turbu- 
lence, but it is not used to derive densities, since the true 
density distribution in observed clouds is unknown. How- 
ever, future measurements of the density distribution can 
be used to apply the 'correction' shown in Figure |4j 



background is the CMB, the 'filling factor' of the molec- 
ular cloud is simply its size relative to the beam size. If 
there is a continuum source in addition to the CMB, the 
size of the continuum source and the intervening molec- 
ular cloud both affect the absorption depth. Throughout 
this paper, we use the term 'filling factor' to refer to the 
fraction of the beam area filled by the absorbing molecu- 
lar cloud and 'covering factor' to refer to the fraction of 
the background continuum source that is covered by the 
intervening molecular material. 

The VLA archival images were used to estimate the 
size of the illuminating background source. When im- 
ages at both wavelengths were available, we separately 
determined the 2 cm and 6 cm source sizes. The source 
size determination is imprecise because we select a sin- 
gle source size for non-uniform surface brightness pro- 
files, and in many cases the VLA observation did not 
recover the full flux densit y seen in single-dish measure- 
ments. Araya et al. (2002) estimated the interferometer- 
to- single-dish flux ratio at 6 cm in this sample and found 
that the interferometer observations recovered anywhere 
from 3% to 100% of the single-dish flux. We repeat these 
measurements at 2 cm and find that the typical recovery 
fraction is higher, ~ 40% to 100%, although sources for 
which only VLA upper limits could be measured have 
recovery fractions < 1%. 




LVG n(H 2 ) cm- 



Figure 4. The mean density from a lognormal density distribu- 
tion plotted against the density derived assuming a single density 
per region (i.e., the directly LVG-derived density). At low densi- 
ties, the wider turbulent distributions are heavily biased towards 
"observing" higher densities than the true mean density. The dis- 
tributions cut off at the low end where the optical depth ratio be- 
comes a double- valued function of density; at these low densities, 
no detections are expected at our survey's sensitivity. The cutoff 
at the high end is where the optical depth ratio becomes constant. 



4. ANALYSIS 



4.1. 



Measuring Line Optical Depth 

In order to measure physical properties of an absorb- 
ing source, measurements must be obtained of the optical 
depths of both the lio — In and 2n — 2 i2 lines. These 
measurements are presented in Tables T] and [I] Once an 
optical depth with errors is determined, the spectral line 
depths can be matched to large velocity gradient (LVG) 
models to determine column and spatial density. The 
spectral line optical depth depends both on the nadir 
flux density of the absorption line and the strength of 
the illuminating background continuum source. If the 
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Figure 5. The filling factor corrected (FFC) density vs. the de- 
rived density with no filling factor correction. While there are some 
cases where the correction results in an order of magnitude or more 
increase in the density, most points show a small correction. The 
black line is the one-one line. Red squares show where the filling 
factor corrected point was used, while blue circles show where the 
uncorrected point was used. Magenta left-pointing triangles are 
limits where the filling factor correction was used, green downward 
triangles are limits where the uncorrected points were used, and 
orange upward triangles are lower limits where the filling-factor 
correction was used. 



The optical depth measurements were "filling factor 
corrected" by assuming the CMB only contributed flux 
density over the same area as the H II region (i.e., the 
foreground cloud covers the exact same patch of sky as 
the UCH II region). When the H II region is small (e.g., 
10% of the beam area or less), the contribution of the 
CMB to the continuum is negligible, but in cases of more 
diffuse H II regions, the CMB contribution is significant, 
particularly at 2 cm. The inferred optical depths and 
source areas are presented in Table [4j Both "filling fac- 
tor corrected" and uncorrected densities are presented in 
Table^ The effect of the filling factor correction (FFC) 
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on density measurements is shown in Figure [5] In a few 
cases, no volume density- column density parameter space 
in the models ( Section [3| was consistent with the spectral 
line ratio after filling factor correction: in these cases, the 
filling factor correction was not used. Similarly, no filling 
factor correction was applied to sources without detected 
continuum. These exceptions are noted in Table [5] in the 
"Flag" column. 

The above definitions are summarized briefly in the 
following equations: 

Sv,obs = S UyCon t(l — CFe Tu ) — S Vi cmb{FFc Tu ) 

FF = £1 cloud I ^6 earn 

C F = Qcloud/ ^continuum 



in which CF is the "covering factor", FF is the "filling 
factor", and there is no positive contribution from the 
CMB because it is assumed to be removed by position- 
switching. 

The systematic uncertainties in the continuum and the 
filling factor result in similar errors in the optical depth 
measurement, and together dominate the total error bud- 
get for our measurements. A 30% error in the lio — In 
and 20% error in the 2n — 2i2 continuum levels were 
assumed because of flux calibration uncertainty charac- 
teristic of the instruments. An additional 10% error in 
the beam area, which sets the maximum coupling to the 
CMB (assuming a beam-filling source), was included to 
account for focus error. A 20% statistical error in the 
cloud filling factor was assumed for the majority of the 
survey, but it was decreased to 10% when the ratio of 
continuum to CMB flux was > 0.5 and the source size 
was small, indicating that the VLA-measured source is 
indeed the dominant continuum component in the beam. 
The statistical error does not account for systematic er- 
rors in the geometric assumptions. Note that changes 
to the filling factor should have a minimal effect on the 
derived density unless the source sizes at 2 cm and 6 cm 
differ substantially, while changes in the filling factor will 
always have a large effect on the derived column density 
(Figure [6]). 

Measurements of volume and column density were 
taken by averaging over the regions of LVG model pa- 
rameter space consistent with both spectral line optical 
depth measurements to within la. The "la" (68% confi- 
dence; errors are non-gaussian) error bars on the derived 
parameters (TV, n, X) were taken to be the extrema of 
these regions. An example of this fitting process is shown 
in Figure [7] A second example demonstrating a lower- 
limit on the density (instead of a direct measurement) is 
shown in Figure [8] This method is not as robust as y 2 
fitting, but because there are no free fit parameters, a 
statistically meaningful xl cannot be computed. 

In some cases, the ratio of the spectral line op- 
tical depths was consistent with low density (n < 
100 cm" 3 ) and high abundances (X(o-H 2 CO) > lO" 8 
per km s _1 pc _1 ), but these were ruled out based on the 
prior assumption that extremely high H 2 CO abundances 
should not be observed at very low densities, since it is 
formed at higher densities and destroyed by hard UV a t 
low columns (see discussion in Troscompt et al. 2009a). 




Figure 6. The dependence of derived parameters on the filling fac- 
tor, assuming an optical depth ratio ti 10 _i 1;l /t2 11 -2i 2 =1 (solid), 
2 (dash-dot), or 4 (dashed). The X-axis is the "real" optical depth, 
T\ — i(real) = T\ — i(observed) / FF . Assuming the same filling fac- 
tor correction is applied to both the lio — In and 2n — 2i2 lines, 
filling factor correction will only move the measurements along the 
X-axis of these plots. A decrease in the filling factor requires an 
increase in the true optical depth to maintain a constant appar- 
ent r (observed), which in turn drives up the derived abundance 
and column density while leaving the volume density unchanged 
(except at high optical depths, r > 0.2). 



G32.80+0.19 filling factor corrected 




log(n(H,)) 



Figure 7. An example of the column density - density parame- 
ter space available given measured lio — In and 2n — 2i2 opti- 
cal depths. The dashed lines show abundances log 10 (X(o-H2CO)) 
per km s — 1 pc — 1 . The contours show the regions allowed by the 
measurements of optical depth (lio — In: black, 2\\ — 2\2' grey, 
ratio: dotted); the middle curve is the measured value, while the 
pair of curves around it are ±1<7 including systematic error. The 
shaded region shows the allowed parameter space from which the 
physical parameters are derived. 



G33. 13-0.09 Q filling factor corrected 




Figure 8. Same description as Figure [7] but for the strongest 
component in G33. 13-0.09. It was only possible to measure lower 
limits on the volume and column density for this line; it is therefore 
assigned flag 8 in Table [5] 
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4.2. Systematic Errors: Absorption Geometry 

There are potential systematic errors associated with 
geometric assumptions, i.e. the filling factor. There 
are four geometric configurations possible; these are out- 
lined in Table Jm The "small source" geometry (3 and 
4) is technically impossible given that the CMB is al- 
ways present in these observations, but it is equivalent 
to the scenario in which the small illuminating compact 
source (UCH II) is much brighter than the CMB in the 
beam. The second column shows the effects of applying 
the 'true' filling factor correction for errors 2 and 4. For 
error type 3, the optical depth will only be overestimated 
if the absorber is "corrected" to be smaller than the back- 
ground source (i.e., if a correction is applied when none 
should have been). 



Figure [6] shows the effects of incorrect geometric as- 
sumptions. Type 1 and 3 errors - i.e. filling factor over- 
corrections - will result in measurements of column and 
abundance that are greater than the real values, while 
type 2 and 4 errors will result in column and abundance 
measurements that are lower than the real values. 

Additionally, it is possible that an observation will in- 
clude a beam-filling, low-density source that will con- 
tribute negligibly in 2n — 2 i2 line absorption but sub- 
stantially in lio — In absorption over most of the beam 
area. This type of error will result in an underestimate 
of the volume density. 

Since these errors are failures of assumptions, they can- 
not be quantified, but Figure [5] shows the effects of cor- 
recting for these errors to the extent possible with the 
available data. 



Table 6 

H2CO Geometric Systematic Errors 



Real Geometry Assumed filling factor = 1 Assumed filling factor < 1 



1. Beam- filling source, beam- filling absorber a T m = h TR tm > tr 

2. Beam- filling source, small absorber tm < tr tm = tr 

3. Small source, beam- filling absorber tm — tr tm >= tr 

4. Small source, smaller absorber tm < tr tm = tr 



tm = measured optical depth 
b tr = real optical depth 

4.3. RRLs 

Radio recombination lines are used to measure the ve- 
locity of the UCH II regions. The recombination lines 
75-77a were independently fitted with gaussians because 
the signal-to-noise in each spectrum with a detection was 
high. Out of our 24 spectra, there were 21 H detections, 
13 He detections, and 12 C detections; Table [7] shows 
the fitted parameters using the 76a lines (75a and 77a 
were also measured but are not reported for brevity). 
For some of the analysis in later sections, we addition- 



ally u s e the d eeper and more careful RRL study by 



Roshi 



et al. (2005), who observed 17 of our sample in the 89- 
92a lines. We atte mpted to measure carbon RRLs in the 



Araya et ah] ( 2002 ) spectra, who only measured hydrogen 
RRLs. We detected one carbon line in G61.48 and tenta- 
tively (~ 2a) detected another three in G32.80, G34.26, 
and G45.45; we report the low-significance detections in 
these sources because of corresponding detections of C75- 
77a. 

We compare the central velocities of the H and C a 
lines to the velocities of the H2CO absorption lines on 
a case-by-case basis in Figures 23]|40 The spectral line 
profiles are used to fit the ob serva tio ns in to the models 
discussed in detail in Sections HT2l and lol 



Table 7 

Measured RRL 76 properties 



Source 
Name 



H 



Peak Center 
H76a a R76a 
(Jy) (kms- 1 ) 



He 



FWHM 

R76a 
(km s- 1 ) 



Peak 
Ke76a 

(Jy) 



Center 
He76a 
(km s- 1 ) 



FWHM 
He76a 
(km s- 1 ) 



Peak 
C76a 

(Jy) 



c 



Center 
C76a 
(km s- 1 ) 



FWHM 

C76a 
(km s- 1 ) 



G32.80+0.19 


0.622 


15.69 


12.09 


0.066 


16.49 


9.25 


0.015 


15.40 


8.27 




(0.001) 


(0.03) 


(0.03) 


(0.002) 


(0.36) 


(0.38) 


(0.002) 


(1.45) 


(1.64) 


G33.13-0.09 


0.067 


73.49 


14.10 
















(0.001) 


(0.17) 


(0.17) 














G33.92+0.11 


0.157 


101.86 


12.16 


0.013 


99.07 


13.60 










(0.001) 


(0.07) 


(0.07) 


(0.001) 


(0.87) 


(0.87) 








G34.26+0.15 


0.367 


54.68 


10.43 


0.034 


51.98 


6.54 


0.026 


59.54 


5.66 




(0.004) 


(0.06) 


(0.09) 


(0.002) 


(0.46) 


(0.49) 


(0.002) 


(0.55) 


(0.56) 




0.251 


37.46 


22.76 
















(0.003) 


(0.29) 


(0.12) 














G35.20-l.74 


1.016 


47.94 


10.70 


0.105 


48.26 


8.27 


0.045 


44.18 


4.05 




(0.002) 


(0.02) 


(0.02) 


(0.002) 


(0.21) 


(0.21) 


(0.003) 


(0.33) 


(0.33) 


G35.57-0.03 


0.036 


52.38 


13.71 
















(0.001) 


(0.41) 


(0.41) 














G35.58+0.07 


0.044 


46.68 


10.55 


0.007 


43.15 


6.30 










(0.001) 


(0.20) 


(0.20) 


(0.001) 


(0.94) 


(0.94) 








G37.87-0.40 


0.446 


59.99 


15.47 


0.042 


60.16 


11.88 


0.018 


59.27 


7.93 
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G43.89-0.78 


0.103 


54.98 


10.83 


0.010 


54.18 


7.72 


0.007 


54.08 


0.82 




(0.001) 


(0.08) 


(0.08) 


(0.001) 


(0.68) 


(0.68) 


(0.002) 


(0.34) 


(0.30) 


G45.07+0.13 
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~ 
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~ 
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G48.61+0.02 


0.076 


16.77 


10.53 


0.007 


16.33 


7.73 


0.006 
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4.71 




( c\ nni "\ 
(U.UU1 ) 
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G60.88-0.13 


0.067 


18.30 


9.08 








0.023 


21.77 


2.54 




(0.001) 


(0 1 9) 
















G69.54-0.98 


0.017 


3.69 


16.24 
















(0.001) 


(0.64) 


(0.64) 














G70.33+1.59 


0.343 


-19.18 


12.59 


0.032 


-20.14 


10.12 


0.025 


-21.67 


3.28 




(0.001) 


(0.05) 


(0.05) 


(0.001) 


(0.44) 


(0.46) 


(0.002) 


(0.33) 


(0.33) 


G70.29+1.60 


0.545 


-26.97 


17.82 


0.042 


-26.32 


14.53 


0.032 


-24.78 


4.71 




(0.001) 


(0.12) 


(0.09) 


(0.001) 


(0.64) 


(0.68) 


(0.002) 


(0.36) 


(0.42) 




0.066 


-64.41 


13.12 
















(0.002) 


(0.70) 


(0.50) 














G61.48+0.09 


0.566 


25.96 


11.16 


0.046 


28.80 


7.86 


0.059 


21.27 


2.48 




(0.001) 


(0.02) 


(0.02) 


(0.001) 


(0.24) 


(0.24) 


(0.002) 


(0.11) 


(0.11) 



a Some H lines were fit with two gaussian components, in which 
case the second fit component is on the second line below. Errors 
(lcr) are indicated by the numbers in parentheses on the line below 
the measurement. 

5. RESULTS 

5.1. Derived Properties 

The average properties of the spectral line components 
associated with the UCH II regions and the other spectral 



lines representing molecular clouds are shown in Table [8j 
The table includes the mean and median only of spectral 
lines with both lio — In and 2n — 2 12 detections that 
yielded measurements of density; upper and lower limits 
are not included. The full results are presented in Table 

El 



Table 8 

Inferred properties 



UCH II Other Lines (GMC) 

Parameter Median a Mean a RMS a Median b Mean b RMS b KS PTE 



log(H 2 Density) (cm -3 ) 4.95 4.91 0.27 4.49 4.61 0.32 0.022 

log(o-H 2 CO Column) (cm" 2 ) 12.59 12.59 0.44 11.86 11.83 0.20 6x10- 

J(o-H 2 CO) -10.84 -10.80 0.46 -11.16 -11.26 0.45 0.028 



a Spectral line components associated with UCH II regions 

k Other spectral lines (associated with line-of-sight molecular 

clouds) 

There is statistical evidence that the deepest spectral 
line components have higher H2CO column and/or abun- 
dance than the other (GMC) components (Table [8]). It 
is unlikely that this difference could be caused by under- 
estimates of the optical depths in the GMC components 
(type 2 and 4 errors, see Table [6| because the filling 
factor correction should tend to cancel out these errors. 



However, it is possible that, in those cases where the H II 
emission and the CMB emission in the beam are the same 
order of magnitude, type 1 errors have occurred: the H II 
region absorber is much larger than the H II region and a 
significant fraction of the spectral line depth comes from 
absorption against the CMB; this error should have little 
effect on the derived density (see Figure [6]) but may lead 
to overestimates of the derived column density. 

Each identified Gaussian component was associated 
with an UCH II region if it was within 5 km s _1 of the 



H 2 CO toward UCH lis 



15 



RRL peak, since RRLs are assumed to be generated in 
the UCH II regions. Any spectral lines blended with the 
UCH II H2CO lines were also associated with the UCH II 
region. Other velocity components, including those with- 
out corresponding RRL detections, were assumed to be 
from GMCs along the line of sight or part of the larger 
cloud not directly associated with the UCH II region; 29 
components were associated with UCH II regions and 46 
were associated with unrelated line-of-sight GMCs (Ta- 
blej3|. 

The density difference between the two populations 
is significant by a Kolmogorov-Smirnov (KS) test with 
~ 2% probability of being drawn from the same distribu- 
tion (the 'probability to exceed' or PTE in Table |gj). This 
result is in contradiction to the results of Wadiai et al.l 



( 1988| ), who found no significant density difference be- 



tween "warm clouds" and "cold clouds" selected and ob- 
served in the same manner (though with larger beams). 
The difference is likely because the larger beam sizes in 
their study and a failure to include the continuum contri- 
bution of the CMB (which is more substantial in a larger 
beam, especially at 2 cm) , resulting in a type 3 error and 
an underestimate of density for their "warm clouds" in 
particular. 

The measured H2 densities do not display any trend 
with heliocentric distance over the range 2-14 kpc, con- 
trasing with mm-continuum surveys of star forming re- 
gions that tend to measur e lower densities at greater 
distances ( |Reid et al. 2010). The lack of correlation in 
Figure [9] demonstrates the strength of the H2CO densit- 
ometry method: the properties of star-forming gas can 
be explored throughout the galaxy with distance bias 
largely removed. Similarly, no trend with Galactocentric 
distance was readily apparent. 



6.0 
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Figure 9. Derived density plotted against kinematic distance. No 
trend is obvious, demonstrating that the H2CO densitometer is not 
biased by source distance. Black squares represent GMCs along the 
line of sight; red triangles represent UCH II regions. 



Densities were measured within a range 10 4 cm 



-3 < 



n(H2) < 10 6 cm -3 due to sensitivity cutoffs at low den- 
sities and thermalization of the sp ectr al line ratio (ratio 
—> 1) at high densities (see Section 6.5 for a discussion of 
the limitations of the densitometer). On the high density 
end, a lower limit on the density remains interesting, as 
densities n(H.2) > 10 6 cm -3 are close to those of low-mass 
protostellar cores and are a strong indication of runaway 
gravitational collapse, since such high densities are rarely 
observed in non-star- forming regions. On the low density 



end, it should be possible to detect the 2n— 2 12 transition 
with sensitivity improvements ~ 2 — 10 x, a consideration 
that will govern the allocated time-on-source for future 
2n — 2 12 observations. 

5.2. Free-free Contribution to 1.1 mm Flux Density 
Measurements 

It is expected that all young star-forming regions 
should be dust-rich and therefore bright at 1.1 mm. We 
therefore compare the BGPS 1.1 mm, GBT 2 cm, and 
Arecibo 6 cm continuum measurements for sources cov- 
ered by the BGPS in Figure [lOl For a flat-spectrum 
(a « -0.1, Tff « 1; I Wilson et al.||2009| ) free-free con- 
tinuum source, the 2 cm riux density should be 1.34 x 
the 1.1 mm flux density. For an optically thick source, 

5 1.1mm ~ 330 S^cm- 

The objects targeted in our survey include 9 of the 
13 brightest (S. , > 1.5 Jy) sources in the range 

° v l.lrara,40 J ' 

32 < i < 48, and 11 of 26 with S, , W / > 1.0 Jy. We 

' l.lrara,40 J 

use flux density measurements from the 40" apertures 
in the BGPS catalog because they are most appropri- 
ate for determining peak brightness of point-like sources 
(Rosolowsky et al. 2010[ ). Out of the sample within the 
BGPS survey area, 6 of 15 sources have free- free frac- 
tions of at least 30%, but potentially much higher if the 
free-free emission is not optically thin. Since the sam- 
ple was selected from well-known UCH II regions, these 
(rather incomplete) statistics are a warning that most 
of the brightest 1.1 mm emission sources in the BGPS 
are likely to be active UCH II regions and therefore may 
include a significant contribution from free-free emission 
to their measured flux densities (Figure 11). The same 
warning applies to other mm-wavelength galactic plane 
surveys, though the contamination should be less severe 
at shorter wavelengths. 

In order to evaluate the impact of this conclusion on 
the BGPS, we examine the flux distribution of 6 cm con- 
tinuum sources from the MAGPIS survey compared to 
the BPGS in the same area, 5 < i < 42 and \b\ < 0.42, 
which is the full range of the MAGPIS survey excluding 
the galactic center, where the BGPS catalog follows a 
different flux distribution ( |Bally et al.||2010| ). 

In Figure [l2j we plot histograms oTTFe MAGPIS 6 
cm flux density and the BGPS 40" aperture flux den- 
sity along with the best-fit power-law distribution linep^] 
Since the 6 cm power-law distribution is shallower than 
the 1.1 mm distribution, the 6 cm sources can dominate 
at high flux densities, although the power-law fit for the 

6 cm sources significantly overpredicts the highest-flux 
bins and therefore the power- law is not an acceptable 
fit above Se cm > 1 JyF 3 ] The dashed line in Figure 12 



12 The pow e r law was fit using the python translation of the 
|Clauset et al.|(2009) power-law fitter provided at http: //code . 
google . com/p/agpy/wiki/PowerLaw The fitter computes the max- 
imum likelihood value of the power-law a and the cutoff of the dis- 
tribution, below which a power law is no longer valid either because 
of incompleteness or a change in the underlying distribution. 

13 We have tested the consistency of the two data sets with a low- 
c utoff power-law dist ribution by the Monte-Carlo process described 
in |Clauset et al.| ( |2009] ) . The BGPS 40" aperture flux densities are 
consistent with a power-law distribution at the p = 0.64 level, while 
the MAGPIS 6 cm fluxes are inconsistent, with p < 0.001 (where p 
measures the probability that the data are drawn from a low-cutoff 
power-law distribution) 
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Figure 10. Bolocam 1.1 millimeter flux density versus the cm continuum flux density at 2 cm (left) and 6 cm (right). The BGPS 1.1 mm 
flux density is moderately correlated with both cm continuum measurements; the legend shows the regression parameter. The expectation 
for optically-thin free-free emission ( a = —0.1, dotted) and for intermediate spectral index emission (a > 0, dashed) are shown to illustrate 
that some sources have significant free- free contributions at 1.1 mm (the optically thick case is not shown for either 2 or 6 cm because it 
does not fit on the plot). The legend shows the predicted flux densities for a given spectral index en, the regression parameter r, and its 
likelihood p. The brighter sources are likely to be less optically thick in the free-free continuum than the faint sources. 

5.3. Distances 

We measure a kinem atic distance to each source using 
the Reid et al. (2009) rotation curve. We resolved the 
Kinematic Distance Ambiguity (KDA) towards each line 
of sight us ing a variety of met hods described below. The 
method in |Sewilo et al.| ( |2004| ) allows a resolution in favor 
of the far distance lor LJCH II regions with an interven- 
ing molecular absorption line at more positive velocities 
in the first Galactic quadrant. Associations with infrared 
dark clouds (IRDCs) can resolve the KDA in favor of the 
near distance. We compare our KDA resolutions to An- 
|derson et al.| ( |2009[ ) , with whom we agree on all common 
sources except lor G33. 13-0.09, which we place at the far 

method. The 
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Figure 11. The distribution of free- free contributions to the 1.1 
mm flux density assuming the UCH II region is optically thin at 
2 cm, fff = (*S , 2cm/l-34)/»S'i.i mm . While 9 sources are either 
dust-dominated or optically thick at 2 cm, 6 sources have free-free 
contributrions of 30% or greater. The other sources in the sample 
were missing 1.1 mm flux density measurements because they are 
outside the BGPS survey area. 



shows the best-fit power-law distribution of the MAGPIS 
flux densities scaled down by 0.67, which is the expected 
decrement for an optically thin free-free source from 6 
cm to 1.1 mm (spectral index a = —0.1). 

Figure [12|q shows the ratio of the BGPS to the MAG- 
PIS best-fit power-law distribution, indicating that the 
free- free contamination fraction is only large (~ 10%) at 
values much greater than the valid range of the 6 cm 
power law fit, which overpredicts the number of sources 
at S^cm « 1 Jy. However, if any of these sources are not 
optically thin at 6 cm, this fraction could be much larger. 
Additionally, these numbers only describe the sources in 
which all of the 1.1 mm flux is free- free emission; the im- 
plication remains that a large number of 1.1 mm sources 
have a substantial (if not dominant) free-free contribu- 
tion. 

Finally, we emphasize that unless a large fraction of 6 
cm sources are optically thick in free-free continuum, the 
lower flux-density BGPS dust-continuum sample should 
be negligibly contaminated by free-free emission sources, 
but the brightest BGPS sources may have a significant 
free-free contribution. 




i except : 

distance based on the |Sewilo et al. 
derived distances are listed in Table 

5.3.1. Size Estimates 

We estimate the source size using two methods. First, 
we use the VLA mea sure ments of UCH II region sizes. 
As stated in Section |4.1[ the VLA size measurements 
are very uncertain and are simplifications of an evidently 
complicated geometry. We estimate a spherical radius 
r = yj area/ir. Second, we assume the gas traced by 
H2CO and the BGPS 1.1 mm images are the same and 
get a 'size scale' r = 2A r mm (H2)/n(H2) where n(H.2) is 
derived from the H2CO line ratio. 

The sizes derived from the two methods are plotted 
against each other in Figure 13 The sizes estimated from 
the two different methods are not well correlated and dis- 
agree by around an order of magnitude in most sources. 
The disagreement could be because of poor VLA-based 
size estimates, substantial 1.1 mm emission from low- 
density gas, or incorrect dust temperature or opacity es- 
timates. While additional line-of-sight GMCs could in 
principle contribute to the N/n size estimate, the dis- 
agreement for sources even without associated GMCs 
prevents this hypothesis from fully explaining the dis- 
agreement. Therefore, any quantities derived from the 
size - i.e. mass, which depends on r 3 - are even less 
constrained. We therefore do not derive any quantities 
dependent on the intrinsic source size. 

6. DISCUSSION 
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Figure 12. Left: Histograms of BGPS 1.1 mm 40" aperture flux densities (red) and the MAGPIS 6 cm flux densities (black), and their 
respective best-fit power-law distributions (a(l.lmm) = 2.41 db 0.03, a(6cm) = 1.72 db 0.03). The dashed black line shows the MAGPIS 
best-fit power-law scaled down to the expected flux density at 1.1 mm assuming all sources are optically thin. Both distributions appear to 
be reasonably well- fit by power-laws above a cutoff (presumably set by completeness), although the power-law significantly over-predicts 
the number of sources with SQ crn > Uy. The histograms are binned by 0.1 dex, and while the best-fit a and x m i n values are independent 
of the binning scheme, the normalization is not. Right: The ratio of the number of MAGPIS 6 cm sources to BGPS 1.1 mm sources as 
a function of flux density for the best-fit power laws. Only 10 1.1 mm sources are detected above 5 Jy (in 40" apertures), so even the 
brightest detected 1.1 mm sources are not purely free- free, but they probably have a substantial free-free component. 

molecular gas should be ~ 1. 

It is therefore interesting to note that Arp 220, possibly 
the most extreme nearby example of a starburst galaxy, 
has nearly the same column per channel as the low end 
of the UCH II regions, suggesting that it is analagous to 
a scaled-up UCH II region to within a factor of a few; 
the measured density in Arp 220 is consistent with only 
the highest-density UCH II regions. M82, on the other 
hand, has a bright continuum background analagous to 
an UCH II region, but a correspondingly low filling fac- 
tor, implying that it consists of many compact but bright 
sources with a total filling factor 0.001-0.1. Alterna- 
tively, the density and column measurements are consis- 
tent with M82 being dominated by quiescent GMCs, but 
that is unlikely given the starburst nature of the galaxy. 

The gravitational lens source B02 18+357 is a differ- 
ent scenario. Its low density is consistent with that of 
a non-star- forming GMC, while its column per km s _1 
is comparable to the Galactic sample. This source is 
therefore likely to be a sightline through a 'normal' qui- 
escent molecular cloud in its host galaxy, similar to the 
narro w beam of an UCH II re gion through the Galactic 




io J 

UCHII radius (pc) 

Figu re 13. A plot of the two derived sizes discussed in Section 
|5.3.1| The two size estimates are at best very weakly correlated. 
Because of the substantial disagreement between the two methods, 
we choose not to explore any parameters with a strong dependence 
on the size. The plotted point size indicates the number of as- 
sociated line-of-sight GMCs, which in principle could lead to an 
overestimate of the N/n size because of additional mass included 
in the 1.1 mm continuum measurement. 

6.1. Comparison to extragalactic observations 

We compare our measured column and volu me den- 
sities t o a se lection of starburst galaxies from Mangum 



et al. (2008) in Figure 14 All of the extragalactic ob- 
servations have much lower column densities per km s _1 
than we measure in the main lines of most UCH II re- 
gions, but similar volume densities. This discrepancy 
can be easily explained by a difference in the area fill- 
ing factor of molecular clouds in observations of galaxies 
and UCH II regions. In a galaxy, the total area filling 
factor of molecular clouds per km s _1 (which is similar 
but not identical to the volume filling factor) is likely to 
be < 1, even in extreme starbursts; although the galaxy 
may appear to be uniformly filled with molecular gas 
in projection, at any given velocity it is likely to have 
significant gaps of ionized or neutral atomic gas. In con- 
trast, an UCH II region should be completely embedded 
in a molecular cloud that is much larger than the free- 
free emitting continuum region, so the covering factor of 



disk. Zeiger & Darling ( 2010| ) note that there is a range 
of covering factors cited in the literature, which can af- 
fect the measured density and column, but should not 
affect the conclusion that the B0128+357 cloud's den- 
sity is not consistent with that of massive-star forming 
regio ns. The low- density g as is detected partly because 



the Zeiger fc Darling] ( 2010) data are 3.5 x more sensitive 
than ours with a background continuum source of similar 
brightness. 

6.2. Line Profiles 

Despite the many systematics discussed above that can 
affect H2CO absorption measurements with a compact il- 
lumination source, it is possible to directly compare the 
properties of gas along a given line of sight without most 
of these hindering factors. Since most of our spectra have 
kinematically resolved spectral line profiles, it is possible 
to make many density measurements at different veloci- 
ties towards each source. An example of this type of anal- 
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Figure 14. Comparison of the UCH II sample (blue circles are 
measurements, blue triangles are lower limits on volume density 
with poorly constrained column densities), the GMC sample (red 
squares), secondary lines associated with UCH II regions (black 
stars) and the extragalactic sample of Mangum et al. (2008) (green 
squares). The errorbars on the Galactic data points are excluded 
for clarity. The observed galaxies have similar densities to the 
Galactic UCH II sample, but significantly lower column densities, 
suggesting that the molecular gas in these galaxies has a filling 
factor << 1. The lack of direct density measurements of UCHII 
regions at high densities is due to the presence of a dominant back- 
ground source; in Arp 220 a direct measurement of density was 
possible because H2CO was seen in emission. 



ysis is shown in Figure [T5| 

ofairialysis is shown in Figure [I6j 



An example demonstrating 

the need for this type 
in which two lines well-fit by gaussian profiles nonethe- 
less display a density gradient because the line centers 
are significantly offset; the figure also demonstrates that 
the offset cannot be accounted for by any instrumental 
effects. 

Of our sample, 18 of the 24 observed lines-of-sight had 
high enough signal-to-noise spectra (S/N> 5 in at least 
four adjacent 0.4 km s _1 channels in both lines) to mea- 
sure the density in many velocity bins. Of these, 12 have 
different peak velocities in the lio — In and 2n — 2i2 
lines, indicating density gradients in the molecular gas 
with velocity. Figure pJjfr is an example density- velocity 
plot. 

We have classified each high S/N spectrum as gradi- 
ent, envelope, or single based on spectral line morphol- 
ogy. The gradient classification was used for gaussian or 
nearly gaussian lines in which the lio — In and 2n — 2 12 
line centers were offset, indicating a gradient in the den- 
sity with velocity; the color listed in the table indicates 
the direction of increasing density. The envelope classi- 
fication was used for flat profiles on the wings of deeper 
gaussian lines. The single classification was used for lines 
where the lio — In an d 2n — 2 12 velocities matched. Low 
S/N spectra were not classified. Classifications are given 
in Tabled 

Of the 12 sources with density gradients, 6 show an in- 
creased density towards the red and 5 towards the blue. 
One source, G45.12+0.13, shows a slight increase to- 
wards the red over a broad (8 km s _1 ) velocity range, 
but a sharp increase towards the blue over only 1 km s _1 
and is therefore classified as other. 

Figures [23][40] show the ' main line' (associated with 
the UCH II region) profile and the associated density, 
column, and abundance velocity profiles. The density, 
column, and abundance measured for each main line via 
the gaussian fit technique are shown overplotted on the 



profiles with blue squares. In all cases, the gaussian fit 
measurement of density is consistent with the individual 
channels nearby and the gaussian fit measurements of 
column and abundance are consistent with the peak col- 
umn and abundance. The consistency of adjacent veloc- 
ity bins confirms the validity of associating gaus sian com- 
ponents in observations of whole galaxies (e.g., Mangum 
et al.|2QQ8 ) or kpc-scale regions, since on these scales the 
lio In and 2 11 — 2 i2 lines should be blended by kine- 
matics to have the same shape. 

6.3. Comparison of RRLs and H2CO lines 

We compare the density spectra with the fitted RRL 
centroids and attempt to interpret these spectra in the 
context of various simple models of H II region interac- 
tion with molecular clouds. The simple models described 
below may actually be short-lived but recurring stages in 
the normal life cycle of a collapsing clump that i s forming 
massive (M > 1OM ) stars flPeters et aL^OlO ). 

We consider five simple models ol embedded UCH II 
regions. For each scenario, we include a brief descrip- 
tion of the model and an analysis of the observational 
consequences in terms of C and H RRL velocities and 
velocity-density structure. We assume that the carbon 
RRLs are only detected if seen in the foreground of a 
bright source. This assumption is based on predictions 
that C RRLs will be amplified by an order of magni- 
tude even in the presence of a weak background (Natta 
et al. 1994). It is backed by a strong correlatio n be- 
tween the continuum and the C RRL intensity (Roshi 
let al. 120051). We also assume that lower- frequency RRLs 
will have a stronger stimulate d emission component than 
higher- frequency RRLs ( jLockman fc Brown||1978[ ). All 
H2CO absorption is assumed to be against the UCH II 
region in this section. The scenario that describes a given 
spectrum is listed in the figure caption for each spectrum 
and in Table 03 

• SCENARIO 1: STATIC 
In a uniform medium with no bulk motions (i.e., 
no collapse), a massive star ignites and generates 



an expanding H II region. Figure 17 



1. Lower frequency RRLs are blueshifted 
from higher-frequency RRLs because of an 
i ncreased stimulated em ission component 
dLockman fc Brown|[l978l ) 

2. A carbon RRL should be seen at the same 
velocity as or blueshifted from the hydrogen 
RRL line center. 

3. Molecular gas closest to the H II region should 
have the highest density because of compres- 
sion by the expanding H II region. It will be 
at a similar velocity or blueshifted from the H 
RRLs. 

• SCENARIO 2: COLLAPSE 
A massive star ignites while spherically accreting 
from a molecular cloud undergoing bulk (inside- 



out) collapse. Figure 18 



1. The H 2 CO-measured density should peak at 
the velocities most redshifted relative to the 
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Figure 15. Plot of the derived parameters per velocity bin for the main line of G32.80+0.19; the full spectrum is shown in Figure m The 
density peak around 16 km s _1 is slightly redshifted of the H and C RRL velocity centers, although the C RRLs are blueshifted ofthe H 
RRLs, indicating that the PDR has been accelerated towards us along the line of sight. The b lues hifted emission tail is suggestive of an 
outflow. This source cannot therefore be easily classified under any of the scenarios in Section |6.3| but is consistent with components of 
scenarios 2 and 3. 

a. The spectra of G32.80+0.19. The GBT 2n — 2i2 spectrum (red solid) has been smoothed to a resolution of 0.38 km s _1 to match the 
Arecibo (black das hed) spectral resolution . Labe led vertical bars indicate the measured velocity centers of H and C RRLs from this work, 
|Roshi et al"] ( |2005| >, and |Churchwell et al.| ( |2010| ). 

o. The measured densities m each spectral bin. The Y-scale is in logio units. Error bars include a 10% systematic uncertainty in the 
continuum and therefore errors in adjacent channels are not independent. Limits are indicated by triangles. Bins with no information 
above the \-o noise cutoff are left blank. The increase of density towards higher velocities led us to classify this source as a red gradient in 
Table El 

c. The measured column densities per spectral bin. Because these column densities are derived from a large velocity gradient code, they 
are in per km s — 1 pc — 1 units. 

d. The measured abundances per spectral bin. The column and abundance are somewhat degenerate, but it is possible in some cases to 
place tight constraints on the total 0-H2CO column while only placing upper limits on abundance and density. The abundance must also 
be interpreted per km s _1 pc _1 . 

In plots b through the blue square with error bars represents the measured value from Table [5] using gaussian fits to the lines. 
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Figure 16. Comparison of G70.29+1.60 (top) and G70.33+1.59 
(bottom) spectra as observed by Arecibo (black) and GBT 
(red/grey). Note that in G70.29+1.60, the 2n - 2 12 line is shifted 
towards the blue of the lio — In line, while in G70.33+1.59 the 
line centers match well. 




Figure 18. Scenario 2: An UCH II region forms from a gravita- 
tionally unstable cloud undergoing inside-out collapse. See Figure 
|17| for a complete description of the figure. The highest density 
should correspond to the highest-velocity infall, so the 2\\ — 2i2 
line peak should be redshifted of the lio — In line peak. The 
hydrogen recombination line may align with a low-density cloud 
but should be blueshifted of the infalling gas. The carbon RRL 
should be redshifted from the hydrogen RRL and blueshifted from 
the H 2 CO line. 





Figure 17. Scenario 1: An UCH II region forms and begins ex- 
panding spherically in a uniform density gas cloud. A cartoon of 
the geometry seen by the observer is shown on the left side of the 
figure, with arrows indicating expansion and darkness of the gray 
shading indicating relative density. The white region around the 
central star is the ionized UCH II region. On the right side, a car- 
toon of the relative velocity and width of the RRLs and H2CO lines 
is shown. The relative heights of the H2CO lines is representative 
of the observed density; black is lio — In and red is 2n — 2i2- 
The narrow emission line with a ? above it indicates a possible 
blueshifted carbon RRL; its height has no physical meaning. In 
this scenario, the hydrogen recombination and H2CO lines should 
occur at the same velocity, and the H2CO lines should show rel- 
atively low-density (high lio — ln/2n — 2i2 ratio) and modest 
spectral line widths. A blueshifted carbon RRL may form, but is 
not guaranteed. 



Figure 19. Scenario 3: An UCH II region expanding in a uniform 
medium ejects a bipolar outflow. Presumably t he b ipolar outflow 
comes from a disk-accreting source. See Figure^] for a complete 
description of the figure. The outflow (indicated by the cones emit- 
ting from the central source) should have lower column density 
but could have high or low volume density. It will be observed as 
high-velocity blueshifted absorption in a line wing. Carbon recom- 
bination line emitting regions may be destroyed by the outflowing 
material. As in the simple scenario 1, the hydrogen recombination 
line should be at the same velocity as the molecular cloud. 




Figure 20. Scenario 4: An UCH II region expanding in a uniform 
medium sweeps up and accelerates material that undergoes trig- 
gered star formation. Because the highest-density material is the 
swept up material, it should be the most blueshifted. See Figure 
[l7|for a complete description of the figure. The orange and yellow 
circles are meant to indicate triggered star formation. 



H 2 CO toward UCH lis 



21 



? ? ? 




Figure 21. Scenario 5: An UCH II region is seen behind a high- 
density, turbulent gas cloud. The turbulence drives large spectral 
line widths, while the high density makes the lio — In and 2n— 2i2 
line depths very close. The RRL velocity could in principle be at 
any velocity relative to the foreground turbulent cloud. See Figure 
1 1 7| for a complete description of the figure. In this case, the ?'s 
indicate an uncertain velocity for the hydrogen RRLs; a carbon 
RRL is not expected because the H II region is not necessarily 
interacting with the molecular gas. 

hydrogen RRLs. Inside-out collapse dictates 
that the highest densities should be infalling 
at the highest speeds. 

2. The C RRL velocity should be between the 
H 2 CO and H RRL velocity since the PDR will 
be decelerated by radiation and gas pressure 
from the H II region 

3. Since the accreting star should be at approx- 
imately the rest velocity of the cloud, there 
should be little to no gas blueshifted from the 
RRL velocity 

• SCENARIO 3: OUTFLOW 

An accreting massive star generates a massive out- 
flow with a significant component along the line of 
sight. Figure [l9j 

1. Substantial low-column, low- abundance 
per km s _1 gas should be observed at veloci- 
ties blue of the RRL velocities. Densities can 
range from low to high. Covering factors may 
be low. 

2. No carbon RRL is expected from the outflow, 
though if the flow is accelerated by ioniza- 
tion pressure a C RRL should be observed 
blueshifted of the H RRL velocity. 

• SCENARIO 4: SWEEPING 

An expanding H II region pushes on a low-density 
envelope, possibly triggering a new stage of star 
formation as in the "collect and collapse" scenario. 
This scenario is similar to #1 but with either a 
higher-density envelope or with more gas swept up 
(i.e., #4 may represent a more evolved region). 
Figure [20] 

1. The hydrogen RRLs should be red of the 
dense gas and the carbon RRLs. The expand- 
ing H II region should accelerate the dense gas 
blue along the line of sight. 

2. A low-density envelope should persist at the 
same velocity as the H II region 



• SCENARIO 5: FOREGROUND CLUMP 
A high density, highly turbulent or high mass and 
rotating clump of gas is in front of the UCH II re- 
gion or surrounds it. This physical situation may 
exist in all of the above and provides alternate ex- 
planations for any spectral line wings. Figure 21 



1. Moderate density gas from a molecular cloud 
will result in high column but moderate den- 
sity at the center velocity 

2. Wide wings of high density gas will exist both 
blue and redshifted of the highest-column 
point 

6.4. The Filling Factor of Molecular Clouds 

We have measured the density in 19 line-of-sight molec- 
ular clouds in addition to the 18 measurements of den- 
sities around UCH II regions (we only include measure- 
ments, not limits, in these counts). The measured den- 
sity from the H2CO line ratio can be compared to other 
measures of densi ty, e.g. the mean molecular cloud den- 
sity measured by Roman-Duval et al.| (12010) from the 
BU-FCRAO GRS It is clear from Figure [22] that the 
average density in GMCs is typically ~ 2 — 3 orders of 
magnitude lower than densities measured in our sample 
o f line-of-sight GMCs 



Roman-Duval et al. (2010) point out that the mean 
densities they measure are significantly below the critical 
density of 13 CO, n cr = 2.7 x 10 3 cm -3 , indicating that 
they do not resolve the high-density clumps that make 
up the GMCs. Our data indicate that a typical GMC 
consists of n ~ 3 x 10 4 cm -3 gas (the median of our GMC 
subsample excluding upper limits), substantially higher 
than the critical densi ty of 13 CO. Taking the rat io of the 
median density in the |Roman-Duval et al.| ( |2010| ) catalog 
to that in our sample, we derive a volume rilling factor 
of 5 x 10 -3 of dense gas in molecular clouds. 

We measured an additional 20 upper limits towards 
GMCs, all of which are consistent with high densities 
(n(H2) > 10 4 cm -3 ), but could represent a sample of 
lower density (n(H 2 ) ~ 10 3 cm -3 ) gas, in which case our 
'measurement' of the cloud volume filling factor is biased 
to be too low. In order to test for this bias, we need 
to acquire more sensitive observations of the upper-limit 
systems. However, we continue analysis below based on 
the assumption that the cloud filling factor measurement 
is realistic, i.e. assuming that the density upper limit 
measurements have densities consistent with the other 
observed GMCs. 

Can a medium with supersonic turbulence produce 
the same density measurements without having to in- 
voke high-density clumping? Below about n(H2) ~ 10 5 
cm -3 , measurements of density in a turbulent medium 
are biased towards higher densities, i.e. the densities we 
report may be overestimates for GMCs since they have 
a median density n(H 2 ) = 10 4 49 cm -3 . For turbulent 
density PDFs with logarithmic widths Ci n (p)//n(p) ^ 1-5, 
the overestimate is no more than 0.4 dex, and there- 
fore can only bring the fillin g factor up by a factor < 3. 
As discussed in Section 3^ and Figure [4j a high-density 
tail could create a larger discrepancy (~ 0.5 dex). How- 
ever, at the measured densities, these are extreme up- 
per limits on the 'turbulent correction', and therefore 
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Figure 22. Histograms of the GMC and UCH II subsamples from 
our data plott ed along with the G MC-averaged densities from the 
13 CO |Roman-Duval et al.| ( |2010| GRS measurements arbitrarily 
scaled to tit on this plot. 1 ne measured densities in UCH II regions 
are significantly (by a KS test) higher than densities in GMCs. 
The H2CO-measured densities in GMCs are 2-4 orders of mag- 
nitude higher than volume- averaged densities of GMCs from the 
GRS, suggesting that GMCs consist of very low volume-filling fac- 
tor (~ 5 x 1 0~ 3 ) high-density (n(H 2 ) ~ 3 x 10 4 cm -3 ) clumps. In 
Section [6^4] we argue that the observed difference is most likely not 
a selection effect imposed by the different gas tracers. The GMC 
upper limits shown are 3 — a upper limits, and all are consistent 
with the measured GMC densities. 

(gravo) turbulence alone cannot account for the measured 
densities. 

What clumping properties are required to reproduce 
the observed density? As long as the clumps are all opti- 
cally thin in the H2CO absorption lines, the spectral line 
optical depths and ratio are independent of clumping. 
However, a large number of low-density (n(H2) ~ 10 3,5 
cm -3 ) clumps optically thick in the lio — In line and 
thin in the 2n — 2 12 line would appear to have a higher 
density. This phenomenon could only occur at densities 
< 10 4,5 cm -3 , where the lio — In absorption line is much 
stronger than the 2n — 2 12 line, and column densities 
10 14 cm- 2 > 7V(o-H 2 CO) > 10 13 5 cm" 2 per clump (at 
higher columns, both lines are optically thick; at lower 
columns, both lines are optically thin). Assuming a typ- 
ical H2CO abundance (ortho+para) Xh 2 co = 10 -9 , the 
required spherical clump radius would be ~ 0.3 pc, which 
would be Jeans-unstable at the assumed density and tem- 
perature (40 K) and is therefore unlikely to persist for 
long time periods P] We therefore regard a collection 
of optically thick cmmps in the H2CO lio — In li ne to 
be unlikely; clumps optically thick in both lines are even 
less likely following the same line of reasoning. 

The combination of the observed large spatial scales 
(and therefore low volume- averaged density) of GMCs 
and the high densities measured along essentially arbi- 
trary sightlines through these GMCs suggests that GMCs 
are not consistent with a purely turbulent medium with a 
lognormal density distribution. The observations also re- 
quire a more substantial high-density tail than typically 
seen in gravoturbulent simulations, i.e. they require a 
clumpier medium. 

Alternatively, it is possible that H2CO is chemically 

14 However, the lifetime of such clumps in a turbulent medium in 
which small-scall turbulence supports the clumps against collapse 
is unconstrained. 



enriched in high-density pockets within a turbulent 
medium, which would imply that H2CO observations 
probe different gas than CO. No such mechanism has 
been proposed on theoretical grounds, and the timescales 
for enhancement would have to be very short (inter- 
mittent density enhancement occurs on ti mescales much 
shorter than the dynamical timescale; Krits uk et al.| 
2007), so we regard this possibility as unlikely but in- 
clude it for completeness. 

Another alternative is that the 13 CO systematically 
underestimates the mass or overestimates the volume of 
the cloud, resulting in an underestimate of the cloud den- 
sity. Sub-thermal excitation of 13 CO in the low-density 
parts of the cloud can lead to a n underestimate of the 
mass (Roman-Duval et al. 2010, Section 9.3). Since the 



cloud sizes were derived using an assumed spherical sym- 
metry, but molecular gas is typica lly observed in fila- 
menta ry structures, the densities in Roman-D uval et al.| 
(2010) are likely to be lower limits on the mean density 
in the molecular gas. While both of these factors bring 
the H 2 CO and 13 CO densities into closer agreement, it 
is difficult to quantify these effects. 

6.5. Strengths and Weaknesses of the H2CO K-doublet 
Densitometer 

The dynamic range of a spectral line as a tracer of 
a physical quantity is an important consideration when 
designing an experiment. We have demonstrated only 
a modest dynamic range in density measurements using 
the lio — In and 2n — 2i2 lines in absorption against 
bright background sources: above n(H.2) ~ 10 5 ' 6 cm -3 , 
we are only able to set lower limits on the density be- 
cause the spectral line ratio asymptotes to 1 , and below 
n(H.2) ~ 10 4 cm -3 , the 2n — 2i2 line optical depth drops 
to very low levels (Figure [2| . 

The lower limits on density at n(H.2) > 10 5,6 cm -3 
can only be modestly improved upon by using higher 
K-doublet transitions (e.g. 3i2 — 3i3) when observing 
H2CO in absorption against bright continuum sources. 
However, when observing anomolous absorption against 
the CMB, an additional density diagnostic is the tran- 
sition from absorption to emission at higher densities, 
which expands the sensitivity of the lio — ln/2n — 2i2 
pair to n(H 2 ) « 10 6 5 cm -3 . 

The low-density end can only be probed by more sen- 
sitive observations of the 2n — 2i2 line. Because the 
H 2 CO line depths become negligible below n(H 2 ) ~ 10 3 
cm -3 , the densitometer is not a useful probe below these 
densities. However, at such low densities, even within a 
molecular cloud, it is questionable whether any molecular 
probes are reliable, as even CO will be undera bundant 
in these environments (e.g., [Glover et aL|2010|). 



As noted in Mangum et alTf PQQ8) and Zeiger & Dar 



ling (2010b, the K-doublet densitometer has the advan 
tage that line detection only depends on the brightness 
of the background source and the gas density. It can 
therefore be used nearly independent of distance when 
observing clouds against the CM B or bright illuminat- 
ing background sources. The |Zeiger fc Darling (2010) 
measurements are more sensitive than any presented in 
our study because of longer integration times and the 
selection of a bright illuminating background source de- 
spite their target being at a distance z=0.68. Following 
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this line of reasoning, any bright synchrotron or free-free 
source can be used to sensitively probe the density of a 
line-of-sight molecular cloud in the Galaxy. The observa- 
tion will have angular resolution limited only by the size 
of the background source as long as it is much brighter 
than the CMB in the beam. 

7. CONCLUSIONS 

We have presented a pilot study to measure molecu- 
lar gas densities in clouds along 24 lines of sight in the 
H2CO lio — In and 2n — 2i2 transitions primarily to- 
ward UCH II regions . We have shown that the H2CO 
densitometer is robust within reasonable ranges of tur- 
bulent density distributions, most cloud geometries, and 
different cloud clumping properties. We have presented 
the methodology and discussed the errors intrinsic to the 
H2CO densitometer. 

Gas volume densities were measured toward 14 of the 
24 sources using the best-fit gaussian profiles; density 
limits were measured for the remaning 10. In 18 sources, 
it was possible to estimate the density in each 0.4 km s _1 - 
wide channel centered on the main line. Of these, 12 
showed some sign of a density gradient with velocity, 

5 appeared to have a single-valued density (i.e. only a 
single spectral line component well- fit by a gaussian), 
and one, G69.54-0.98, had a spectral line optical depth 
that was beyond our ability to model. 

Velocity-density gradients have been used to fit 18 
sources with simple models of UCH II regions embed- 
ded in molecular clouds. We have found some examples 
consistent with inside-out collapse onto the UCH II re- 
gion, UCH II region expansion, and bulk outflow. H2CO 
absorption provides a unique probe of the physical con- 
ditions around UCH II regions because it is only seen 
in absorption against the continuum background (for 
sources much brighter than the CMB), giving different 
constraints than mm and sub-mm spectral lines that are 
seen both in front of and behind the UCH II region. 

The measurements of serendipitously detected line-of- 
sight GMCs revealed densities ~ 200 times higher than 
volume-average densities measured using 13 CO. The 
high density measured suggests that GMCs consist of 
many sparsely distributed high-density clumps and have 
density distributions inconsistent with the lognormal dis- 
tribution predicted by supersonic turbulent models. The 
implied density distribution is also more skewed to high 
densities than predicted by typical gravoturbulent sim- 
ulations. Alternatively, the 13 CO-based mean densities 
may be lower than the mean densities within the molec- 
ular gas either because they underestimate the mass or 
overestimate the volume of GMCs. 

The density measurements show that UCH lis are asso- 
ciated with high-density (n(H.2) > 10 4,5 cm -3 ) gas, and 
UCH lis are associated with higher column and volume 
densities than other line-of-sight molecular clouds, in 
contradiction with previous results ( Wadiak et al.|1 988). 

The 6 cm, 2 cm, and 1.1 mm flux density measurements 
are strongly correlated and in most objects in our sam- 
ple the 1.1 mm flux density has a substantial (> 30%) 
contribution free-free continuum emission. This result 
implies that the brightest sources detected in the BGPS 
have significant free-free emission. A comparison to the 

6 cm MAGPIS survey suggests that the sample of 1.1 
mm sources below about 3 Jy is not significantly con- 



taminated by pure free- free emission sources. 

Comparison of the densi ty measurement s in ou r sample 
to the starburst sample of Mangum et al. (2008) suggest 
that the molecular gas volume filling factor m most of 
these galaxies is small (~ 0.01), but in Arp 220 it is 
quite high (> 0.1). The physical properties measured by 
H2CO in Arp 220 are similar to those in UCH II regions. 
Although velocity-density gradients are observed in our 
sample, we argue that kinematic spectral line blending 
should uphold the assumption of a single spectral line 
profile in galaxies as robust for radiative transfer pur- 
poses. 
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Figure 23. Top: The GBT 2n - 2 12 and Arecibo li - In spectra of G33. 13-0.09. Bottom: The GBT H75o; and Arecibo HllOo; spectra 
with the GRS 13 CO spectrum overlaid. The left axis is for the RRLs and the right axis is for the 13 CO. 
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Figure 24. Same panel description as Figure [15] The recombination lines in this source show enormously different velocities. HllOa 
is 14.6 ib 2.9km s" 1 reds hifted from H75-77a, which is a significant fraction of the total line width (<7no = 41.6 db 6.35, (775 = 31.8 db 2.2). 
|Churchwell et al. Il2010l measure H30oi slightly bluer still, 



= 71.4 db 1.7km s -1 ,<73o = 27.8 db 3.2km s -1 . This discrepancy could be 
explained by a significant gradient in the density of the ionized gas with velocity, i.e. low-density gas blowing out from the back side of 
the UCH II region. The tail of redshifted gas from 80-85 km s - ^ is either rapidly-infalling, low-column gas (Scenario 2) or unassociated 
line-of-sight cloud material (Scenario 5). The main line has a number of points in which only a lower limit on density could be measured, 
but because of the rising lower limits out to ~ 78km s — 1 , the spectrum is classified as a red gradient. 
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Figure 25. Left: Same description as Figure [T] Right: Panels described in Figure [15] The red gradient and carbon recombination line at 



systematic error bars included for the main line resulted in a lower limit on density being measured. In this spectrum, a gaussian line fit 
is not representative of the line profile. 



26 



Ginsburg et al 




Figure 26. Le/fc: Same description as Figure^ Right: Panels described in Figure^] The red velocity gradient in this source is very 
steep and widely separated (by 6 — 15km s _1 ) from the RRL velocities, indicating that the dense gas is rapidly approaching the UCH II 
region. The line profile is exemplary of Scenario 2, inside-out collapse, since the dense gas is redshifted of the recombination lines and the 
C RRL is at an intermediate veloc ity. The derived de nsity at the peak is lower than previous measurements of density using C 34 S and the 
mm transitions of H2CO jMangum fc Wootten|1993| >, although both of those measurements are consiste nt with the density in at lea st one 
velocity bin. The measurements are consistent despite the high temperature (Tk = 230K) measured by (Mangum & Woot ten|1993| >. 
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Figure 27. Left: Same description as Figure^ Right: For panel descriptions, see Figure [15] Identified as a single gaussian line. The 
line peak could be optically thick based on morphological comparison to G69. 54-0.98. Since the carbon RRLs are aligned with the H2CO 
column peak and both are blueshifted from the hydrogen RRLs, this spectrum is consistent with Scenario 4, in which a layer of gas is being 
pushed towards the observer and compressed. 
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Figure 28. Left: Same description as Figure^ One of the two off positions contained an absorption line at 45 km s -1 ; the position was 
coincident with an IRDC. Right: For panel descriptions, see Figure [l5| The densities measured are very high (n(H2) > 10 6 cm -3 ) and 
lower limits are all at similar levels. The match between the H2CO and rtRL velocities and the slight blueshift of the lower frequency RRLs 
suggest that this source is consistent with Scenario 1, though the lack of a carbon RRL detection suggests that the UCH II and molecular 
gas are not interacting. 
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Figure 29. Left: Same description as Figure ^ Right: For panel descriptions, see Figure ^] Identified as a blue gradient because 
of the low density region centered around 53 km s — 1 . However, in this source, the low density region is also substantially lower column 
per km s — 1 pc — 1 , so it may be more appropriate to classify as an envelope. Because the peak density is slightly redshifted from the RRL 
velocities, this source is consistent with Scenario 2, but the lack of carbon RRLs and the lower limits on density at many velocities makes 
the match imperfect. 
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Figure 30. Left: Same description as Figure^ Right: Plot of the derived parameters per velocity bin for G37.87-0.40. This source is 
an example with high column and density associated with the line core (55-63 km s — 1 ) and declining columns at the edge of the line core 
(50-55, 63-67 km s _1 ) with additional low density (n < 10 4,8 cm -3 ) components along the line of sight. A modest increase in density 
towards bluer velocities is observed, but the complexity of the line morphology suggests that this gradient may represent unassociated line 
of sight clouds. The relative velocities of the carbon and hydrogen recombination lines are consistent with the very simple Scenario 1, 
although the H2CO line suggests a more complex situation. 
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Figure 32. Left: Sa me d escription as Figure ^ Right: Plot of the derived parameters per velocity bi] 
descriptions, see Figure Il5] The line core (55 km s _1 ) is very deep in the lio — In line but the lio — ln/2n 



bin for G43.89-0.78. For panel 
2i2 ratio, and therefore 

the spatial density, is substantially higher in the blue tail (48-52 km s _1 ). Since the column densities are lower in the blue tail, it probably 
represents a relatively small patch of highly turbulent, high-density material as in Scenario 5. It is also possible that the blue tail represents 
an outflow as in Scenario 3, but the high density measured is not expected in an outflow. It is possible that there is a signficant contribution 
to the 55 km s _1 lio — 111 line from surrounding low-density gas in the beam absorbing against the CMB, creating a Type 3 error (Table 

E 
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Figure 3 3. Left: Same panel description as FigureJD The H75oi line cannot be fit by a single gaussian. Right: For panel descriptions, 
see Figure [15] The velocity difference between the RKLs and the H2CO lines is consistent with Scenario 2, but no carbon RRL is observed 
so there is no direct evidence of interaction between the UCH II region and the molecular gas. 
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Figure 34. Left: Same description as Figure ^ Right: For panel descriptions, see Figure [15] The line is classified as other because of the 
very wide, flat line. The relative velocities of the recombination lines and their alignment with the H2CO line are consistent with Scenario 
1. 
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Figure 35. Left: Same description as Figure [T] Right: For panel descriptions, see Figure [15] Identified as a blue gradient because density 
clearly increases towards the blue, while only upper limits are measured towards the red. Based on the relative locations of the RRLs and 
the peak H2CO column, the spectrum appears to be consistent with Scenario 2. 
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Figure 36. Left: Same description as Figure pi Right: Panels described in Figure ^] Identified as a red gradient, 
recombination line velocities are consistent with [Scenario 1. However, the highest density gas is the most redshifted, 




The hydrogen 
which is more 



consistent with Scenario 2. Additionally, there is a low-density blueshifted tail that is suggestive of an outflow as in Scenario 3. 
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Figure 38. Same description as FigurefT] There were not enough adjacent high S/N velocity bins to perform the velocity-density analysis. 
However, the velocity coincidence of the aetected H RRLs and H2CO lines suggests that this source is an example of Scenario 1, but no 
carbon RRL is detected. 
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Figure 39. Left: Same description as FigurefT] Right: Panels described in Figure [15] It was only possible to derive upper limits on 
density at most velocities. The locations of the RRLs and H2CO lines are consistent with Scenario 2. 
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Figure 40. Left: Same description as Figure [T] Right: Panels described in Figure [15] Identified as a single gaussian morphology because 
of the peak in density and column. Because both the H2CO and carbon RRLs are blueshifted from the hydrogen RRLs, this source is an 
excellent example of Scenario 4. 
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Figure 41. Left: 



Panels described in Figure [l5| 
-i ; c ^rV. ere 



The flat-bottomed line profile suggests that the 



Same description as Figure RH Right: 
absorption in this source is optically thick. The blank region around 8 km s — 1 is where there is no allowed parameter space in the models. 
The rest of the measurements are likely to be incorrect because they are biased by optically thick components contributing to the line. 
Because of the added uncertainty in the density measurement from the high line optical depth, it is not clear what scenario, if any, this 
source's line profile is consistent with. However, the RRL at 16.4 km s — 1 (off the right side of the spectrum) would be consistent with 
Scenarios 4 and 5. 
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Figure 42. Left: Sa me d escription as Figure ^ Right: Plot of the derived parameters per velocity bin for G70.29+1.60. For panel 
descriptions, see Figure^] The velocity offset between the lio — In and 2n — 2i2 line centers results in a strong blue density gradient. 
Since the density peaks at the same velocity as the UCH II region, but there is high column gas redshifted of the density peak, the spectrum 
is not directly consistent with Scenario 2, to which it is otherwise similar. 
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Figure 43. Left: Same description as Figure 1^ Right: Plot of the derived parameters per velocity bin for G70.33+1.59. For panel 
descriptions, see Figure^] Unlike G70. 29+1. 607 there is no offset between the lio — In and 2n — 2i2 line centers, so this source is 
consistent with a single gas component along the line of sight. Because carbon recombination lines are detected, though, it is evident that 
the UCH II region is interacting with the molecular gas. The recombination line velocities are consistent with Scenario 1, but the lack of 
a gradient in the H^CO-derived density is not. The region could be intermediate between Scenario 2 (infall) and Scenario 1 (expansion), 
causing the observed density to be more evenly spread in velocity even though it varies spatially. 
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Left: Same description as Figure^ Right: Panels desc ribed in Figure [15| There is no continuum emission or RRL detected 
source: the absorption is exclusively against the CMB. Urquhart et al.|(|2009|) re 



Figure 44. 

towards this source: the absorption is exclusively" against the CMB. jUrquhart et al.| l^D09) report a HCH II detection at 0.23 mjy at 6cm, 
which is substantial ly lower than the ~ 20 mJy absorption detected m the lio — In line. The presence of a faint HCH II region and water 
masers jSunada et al.||2007| > suggest that this is an active but very young region, probably the youngest in our sample. It is therefore a 
good example ot what can be expected of a survey selected from millimeter sources instead of UCH II regions. Because of the lack of an 
RRL detection and the low S/N, we cannot classify this source in an UCH II scenario. 
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Figure 45. Same description as Figure [l] The 2u — 2 12 noise levels are too great to measure physically interesting limits on the density. 
No continuum emission is detected at 2 cm, 6 cm, or 1.1 mm. There is also no evidence of an IRDC in the Spitzer 8/xm image. This 
sight line is probably more accurately characterized as a lack of emission region rather than a dust-absorbed IRDC. The three lio — In 
absorption features detected are associated with 13 CO clouds and only have density upper limits from our data. These observations are a 
strong indication that deep enough H2CO observations in the Galactic plane will always detect H2CO against the CMB or diffuse galactic 
background, and therefore it can be used to measure (or constrain) densities of even diffuse GMCs. 
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Figure 46. Same description as Figure^ The 2n — 2 12 noise levels are too great to measure physically interesting limits on the density. 



